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1.0 INTRODUCTION 


An elastic-plastic fracture mechanics (EPFM) software module for inclusion in the NASGRO 
computer program has been developed by Southwest Research Institute under NASA Marshall 
Space Flight Center Contract (MSFC) NAS8-37828. These modules will hereafter be referred to 
as the Phase 1 development. The fracture and fatigue assessment capabilities developed in 
Phase 1 and the theoretical basis of the EPFM approach based on the J-integral, are described in 
the final report for that contract (R.C. McClung, G. G. Chell, Y.-D. Lee, D. A. Russell, and G. E. 
Orient., “ Development of a Practical Methodology for Elastic-Plastic and Fully Plastic Fatigue 
Crack Growth ”, August 1998). In particular, the User Manual written in support of the EPFM 
modules is contained in Appendix K of that report. The reader is referred to Appendix K for 
further details regarding the J formulations employed in the calculations, the fracture analyses 
that can be performed, how to run the modules and examples of the input data needed, and 
validation of the code. 

Since the release of the NASGRO EPFM module in 1998, two further enhancements have 
been made, herein called Phase 2 and Phase 3. In Phase 2, performed under MSFC Contract 
H-33940D, “ Practical Analytical Tools for Nonlinear Fatigue Crack Growth ,” and completed in 
March 2002, the library of J-integral solutions was improved. In Phase 3, performed under 
MSFC Contract NAS8-02051, "Proof Test Design and Analysis and completed in September 
2003, software modules for implementing proof test methodologies were developed. 

The Phase 2 enhancements to the NASGRO EPFM modules included the following: 

(1) Extension of the EPFM solutions for surface cracks (Model SC01) to include surface 
cracks in rectangular plates subjected to arbitrary uniaxial stressing (Model SC02) 
and surface cracks on the inside and outside of hollow cylinders subjected to arbitrary 
non-linear hoop stresses (SC04). 

(2) Extension of the EPFM solutions for centrally embedded cracks subjected to uniform 
stressing (Model EC01) to cracks subjected to arbitrary non-linear stresses (Model 
EC02) and improvements in the J solutions for uniform stressing (Models EC01 and 
EC02). 

(3) Improvements in the accuracies of the EPFM comer crack J solutions for bending 
(Model CC01). 

The improved EPFM solutions developed in Phase 2 have been incorporated into the NASGRO 
analysis options 5 (J computations), 6 (failure analysis) and 7 (fatigue life analysis). 

The Phase 3 enhancements involved significant additions to the EPFM analysis capabilities 
of NASGRO to facilitate accurate proof test analyses. The proof test methodologies that 
underpin the proof test modules were developed by SwRI under MSFC Contracts NAS8-37451, 
"A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies, ” and NAS8- 
39380, “ Guidelines for Proof Test Analysis. ” In order to accurately implement these two 
methodologies, the one degree-of-freedom (1-DOF) ductile failure modules for surface cracks, 
comer cracks, and embedded cracks developed in Phase 1 were replaced by two degree-of- 
freedom (2-DOF) modules. 
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The Phase 3 enhancements to the NASGRO EPFM modules are described below: 


(1) Extension of the Phase 1 EPFM ductile failure module in NASGRO from a 1-DOF 
assessment (where the severity of surface, comer and embedded cracks is 
characterized by a single value for the crack-tip driving force, J or K) to a more 
accurate assessment based on 2-DOF (where the severity of surface, comer and 
embedded cracks is characterized by two values for the crack-tip driving force, J or 
K), enabling changes in crack shape during ductile tearing to be more accurately 
modeled. 

(2) Addition of a proof test module to implement the procedures in Guidelines for Proof 
Test Analysis to facilitate the use of these by practicing engineers. This module leads 
the engineer step by step through the various stages needed to perform a proof test 
analysis. The module also incorporates service analysis routines that can be used to 
determine fatigue crack growth lives, and critical crack and critical load routines. 

(3) Addition of a tear-fatigue crack growth module for ductile materials enabling the 
behavior of fatigue cracks growing near instability to be quantified. It is well known 
that near instability, the growth rate of cracks can be greatly accelerated. This routine 
is used in the multiple cycle proof test routine described in (4). 

(4) Addition of multiple cycle proof test (MCPT) reliability analysis module to 
implement the procedures described in A Comparison of Single-Cycle Versus 
Multiple-Cycle Proof Testing Strategies. This module includes a probabilistic 
analysis for taking into account the effect of the distribution in initial crack sizes on 
the reliability of a fleet of components entering service after MCPT. This module can 
be exercised to determine the change in service reliability of a MCPT compared with 
performing no proof test or a single cycle proof test. 

The improved EPFM solutions developed in Phase 3 have been incorporated into NASGRO 
through enhancements to option 6 (2-DOF failure analysis for ductile materials), and the 
additions of options 8 (single cycle proof test analysis), 9 (tear-fatigue analysis for ductile 
materials) and 10 (MCPT reliability analysis). 

A summary of the current capabilities of the EPFM module and the phase under which they 
were developed is provided in Table 1. A schematic of the crack models for which EPFM 
solutions are available is shown in Figure 1. 
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Table 1. NASGRO EPFM options developed in Phases 1, 2, and 3. 


Opt 

Analysis Type 

Phase 1 

Phase 2 

Phase 3 

No. 


Crack Models and Loading 

5 

Elastic-plastic 
J computation 

TC01 - Tension 
TC02- Tension 
Bending 
EC01 - Tension 

CC01 - Tension 
Bending 
SC01 - Tension 
Bending 

TC01 - Tension 
TC02- Tension 
Bending 

EC01 - Tension (Improved) 
EC02 - Arbitrary stress (New) 
CC01 -Tension 

Bending (Improved) 
SC01 - Tension 
Bending 

SC02 - Arbitrary stress (New) 
SC04 - Arbitrary stress (New) 

TC01 - Tension 
TC02- Tension 
Bending 
EC01 - Tension 
EC02 - Arbitrary stress 
CC01 - Tension 
Bending 
SC01 - Tension 
Bending 

SC02 - Arbitrary stress 
SC04 - Arbitrary stress 

6 

Elastic-plastic 
failure analysis 

Critical crack 1-DOF 
Critical load 1-DOF 

Critical crack 1-DOF 
Critical load 1-DOF 

Critical crack 2-DOF 
Critical load 2-DOF 



All Phase 1 Models 

All Phase 2 Models 

All Phase 3 Models 

IN 

Elastic-plastic 
fatigue life 
analysis 

All Phase 1 Models 

All Phase 2 Models 

All Phase 3 Models 

8 

Single cycle proof 
test analysis 



Safe Life Analysis 
Critical flaw size 
Fatigue Life 
Proof Test Analysis 
Proof load 
Flaw screening 
Final crack size 





All Phase 3 Models 

9 

Tear-fatigue 



Ductile Materials 





All Phase 3 Models 

10 

MCPT analysis 



Ductile Materials 





All Phase 3 Models 


3 























This Addendum to the User Manual in Appendix K of the final report Development of a 
Practical Methodology for Elastic-Plastic and Fully Plastic Fatigue Crack Growth (hereafter 
referred to as Appendix K) provides a description of the new analytical developments and 
software modules resulting from Phases 2 and 3, validation of the software modules, and 
examples of applying the new modules. Validation of the developments made under Phase 1 is 
presented in Appendix K. 
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Figure 1. The NASGRO crack models for which EPFM solutions ar? available. 

2.0 TECHNICAL ISSUES RELATED TO PHASE 2 ENHANCEMENTS 
The new analyses performed in support of Phase 2 developments consisted of: 

- extension of the SC02 solutions to arbitrary stressing, 

- addition of EPFM J solutions for axial cracks in cylinders subject to arbitrary hoop 
stressing (SC04), 

- improvement in the accuracy of the J solutions for comer cracks subjected to bending 
(CC01), 

- improvement in the accuracy of the J solutions for embedded cracks subjected to 
uniform stressing (EC01) 

- addition of J solutions for embedded cracks in arbitrary stress fields (EC02). 
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In the J formulation scheme used in the NASGRO EPFM Module, J is expressed as the sum of 
elastic, J e , and plastic, J p , components, J = J t {a tg ,o p + a * )+ J p (a,a p ), where a e g is an 

effective crack size equal to the original size, a, plus a plastic zone correction, a signifies the 
applied stress, and superscripts p and s refer to primary and secondary loads, respectively. The 
enhancements to the J solutions incorporated in NASGRO during Phase 2 were predominantly 
related to the plastic component of J, J p . Therefore, the technical aspects of this Addendum 
addresses specifically those issues related to the computation of J p , as J e (a e ff) can be determined 
using linear elastic fracture mechanics principles. In all cases, the reference stress method 
(RSM), see Appendix K, is used to implement J solutions in the NASGRO EPFM module. In 
this formulation, it is important to remember that J p is a function of the primary component of 
loading, as it corresponds to fully plastic fracture behavior and secondary loads cannot influence 
this behavior. 

In the RSM, the plastic component of J, , is given by an equation of the form 


for tension loading and 


r a a b\ .. 

f P l 



\t c c) 
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t RSM 
J P 


( a a b'] 
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- - - 

fiVa 

- , * 
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( 1 ) 


( 2 ) 


for bending, where V is a dimensionless structural parameter and a, <j 0 » and n are material 
properties defining the Ramberg-Osgood equation describing the uniaxial stress-strain behavior 
(see Appendix K). 

The values of V, P* , and M* 0 can be determined from FEA results for hi using the 

optimization RSM scheme described in “ Development of a Practical Methodology for Elastic- 
Plastic and Fully Plastic Fatigue Crack Growth”. In this optimized RSM approach, values of V 
and an optimum yield load, P* (or an optimum yield moment, M’) are found such that the 

RSM reproduces J values derived from finite element analysis (FEA) as accurately as possible. 
The results of the optimized procedure demonstrate the maximum accuracy that can be obtained 
using the RSM. However, due to the limited number of FEA J solutions that can be generated, it 
is not practical to employ the optimized RSM results directly in the NASGRO J module. 
Instead, a pragmatic approach is followed and the module that uses average values for V and 
approximate equations for P* and M* 0 based on simple plastic limit load analyses and empirical 

fits to the actual derived optimized loads. This pragmatic approach is herein called the hybrid 
RSM. 


Plastic collapse loads are defined as 
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where subscript c signifies collapse and the flow stress is defined as 


m 


(4) 


Failure is predicted when the applied load exceeds the plastic collapse load, irrespective of the 
applied J value. 

3.0 J SOLUTIONS FOR SURFACE CRACKS IN PLATES SUBJECTED TO 
ARBITRARY STRESS FIELDS (Crack Model SC02) 

3.1 Implementation of J Solutions for SC02 

The SC02 J solutions were implemented in the EPFM module via the hybrid RSM method 
(see Appendix K). In order to determine the RSM solutions, existing NASGRO stress intensity 
factor (SIF) solutions for surface cracks in arbitrary stress fields were employed together with net 
section yield loads derived for surface cracks subjected to combined tension and bending loads. 
The latter solutions are needed because, in general, arbitrary stress fields when integrated over 
the load bearing section produce tensile forces and bending moments. 



Figure 2. Schematic of SC02 geometry showing location of the neutral axis under combined 
tension and bending. 

The net section yield load for combined tension and bending for model type SC02, 
characterized by the tensile yield load P* , is derived from a plastic limit analysis assuming a 
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neutral axis midway across the net section thickness. In reference to the cross section defined in 
Figure 2, the variation of net section thickness, t^z), with location z as the plate is traversed is 
given by 




(5) 


In this equation, t is the thickness of the plate, a is the depth of the flaw and 2c its total surface 
length. In the case of pure tension, the net section yield load, Po*, is derived from the load 
redistribution due to area reduction and is given by 


P>J3a c 


bt- 


mc 


( 6 ) 


where W-pb is the total length of the plate and <7 0 the yield stress. For the SC01/SC02 models 
P=2, and for CC01 P=l. In the case of pure bending, the net section yield moment, Af * , can also 
be determined analytically assuming the form of the neutral axis given in equation (5) as 


2 

M* = —(-3atn + 4a 2 +6 t 2 )+-(b-c) 

24 4 


(7) 


In the case of combined tension and bending, a proportionality factor, X, is introduced defined as 




Pt 


( 8 ) 


where M is the applied moment and P is the applied tensile load. In the SC02 model, the values 
for P and M are derived from user specified arbitrary stress distributions by integrating these 
distributions over the area of the plate. 

From plastic limit load theory, the equation for the combined tension and bending yield load, 
p; w , can be written as 


K M=fl 


•v 


AtP 


Kj 


\I/2 


+ 4 


M* 


(9) 


In this equation, P 0 * [X] is a net section tension yield load for combined tension and bend 

loading. It equals the value of the tensile load that causes net section yielding under proportional 
loading in the presence of an applied moment related to the tensile load by the proportionality 
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constant X given by equation (8). The value of Pj (l) reduces to P' 0 and M* in the cases of 
pure tensile loading and pure bending, respectively. Pure tensile loading is defined by X=0, thus 

K = Pi U = 0) = flaiylbi - (10) 


Pure bending load is defined by X=°°, thus 
M]= AtP* (A = °o) = fio L 


2 

— (-3at;r + 4a 2 +6 t 2 )+-(b-c) 
24 ' ' 4 V ’ 


( 11 ) 


In the hybrid RSS method, the plastic component of J, J* SM , for combined loading is given by 
the equation 


jRSM 
J p 


= J. 


a a b 
t c c 


V 1 




A W(A)d 
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c )j 


( 12 ) 


where V(A) is a dimensionless structural parameter for the combined loading. Since the value of 
V(A) was only determined in Phase 1 for the two extreme cases of pure tension, V(A=0), and pure 
bending, V(A= °° ), its value for combined loading is herein interpolated between these two 
extreme values using the equation 

V{A) = V(A = oo)(l - exp" 0693 ^ 1 )+V(A = 0)exp‘ <) - 693W 1 (13) 

i 

The deepest point and surface point values, respectively, of V used in SC01 are 1.0412 and 0.973 
for V(A=oo) and 1.8164 and 1.2561 for V(A-O). These two sets of extreme values for V indicate 
the maximum inaccuracies in V(A) that could be generated using the interpolation equation given 
by equation (13). 

3.2 Validation of J Solutions for SC02 

The SC02 J solutions were implemented for arbitrary stressing by utilizing the existing SC02 
SIF routines, adding a routine for determining the applied force and moment corresponding to 
the arbitrary stress, and introducing a net section yield load solution for combined tension and 
bending. 

The new J solutions were partly validated by applying arbitrary primary loads that simulated 
uniform stressing and bending and comparing the resulting J values with the values obtained 
from running the SC01 model for tension and bending. The results are shown in Figure 3 where 
J values derived from the SC01 model are plotted against J values computed using the new SC02 
model solutions. Perfect agreement between the two sets of solutions occurs when the data 
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points fall on the 1 to 1 line. It can be seen from Figure 3 that excellent agreement obtains 
between the SC02 and SC01 solutions, indicating that the integration routines used to determine 
the external forces and moments from the arbitrary stress distribution specified in SC02 and the 
resulting net section yield solutions are correctly calculated. 

Two additional verification tests for SC02 were performed. In the first, a self-equilibrated 

x ( x^\ 

primary stress of the form 1-6— + 61 — I was applied. This form of stress integrates to zero 

force and moment. The resulting J solutions correctly gave non-zero values for J e and zero 
values for J p . In the second validation exercise, arbitrary stresses were specified that 
corresponded to combined tension and bending and the resulting J values where compared to the 
results of manual calculations performed using a spreadsheet. There was exact agreement 
between the two sets of results (see Appendix 4). 



Figure 3. Comparison of J estimations obtained using SC01 and SC02. The arbitrary stress 
fields used in the SC02 computations were chosen to simulate uniform tension and 
pure bending, respectively. The deepest point corresponds to the a-tip and the surface 
point to the c-tip. 


4.0 J SOLUTIONS FOR AXIAL SURFACE CRACKS IN CYLINDERS SUBJECTED 
TO ARBITRARY HOOP STRESSES (Crack Model SC04) 

4.1 Implementation of J Solutions for SC04 

The SC04 RSM J solutions were implemented in the EPFM module using existing NASGRO 
SC04 SIF solutions for internal and external axial surface cracks. The net section yield pressure, 
P * , was taken as equation (B4.9) in Guidelines for Proof Test Analysis, which is a modified 
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form of the equation given by Keifner, Maxey, Eiber, and Duffey in Failure Stress Levels of 
Flaws in Pressurized Cylinders (ASTM STP 536, pp. 461-481). This equation is 


(14) 


In this equation, R m is the mean radius of the cylinder, /?,• is the inner radius, and 

A/A-f|=fl + 1-05VV ,p = - ~ (15) 

'M t J jRJ 

The RSM for the plastic component of J for the SC04 geometry is 



where P is the internal pressure, derived by integrating the hoop stress through the wall of the 
cylinder. I 

•j 

It is important to note that the primary (internal pressure) load can be input in the EPFM 
module in two ways. In the first, the user specifies the actual pressure and the program internally 
determines the hoop stress distribution corresponding to that pressure. For external cracks, the 
derived hoop stress distribution is used in the SIF calculations. For internal cracks, in order to 
allow for the effects of internal pressure acting on the crack faces, the pressure is added to the 
derived hoop stress and this combined stress field is employed in the SIF calculations. For both 
internal and external cracks, the user specified pressure is used in the evaluation of J p with 
equation (14) used for P * . 

In the second method of defining the applied load, the user directly specifies the hoop stress 
distribution through the wall. In this case, for external cracks, this stress distribution is used to 
determine the SIFs and the integrated stress through the wall to determine the internal pressure 
corresponding to this distribution. For internal cracks, it is assumed that the user specified hoop 
stress includes a uniform stress component equal to the internal pressure . This stress distribution 
is used in the SIF calculations to allow for the effects of the internal pressure acting on the crack 
faces. However, the effects of pressure on the crack faces is not included in the determination of 
J pt and equation (14) is used for P* with the pressure P evaluated by integrating the user 
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specified hoop stress and multiplying the resulting “pressure” by the factor R/R 0 in order to 
obtain the actual pressure, P. 

FEA J data is not available for pressurized pipes to allow evaluation of V, so, in lieu of more 
accurate values, V for the deepest and surface crack positions are both set to 1. 

4.2 Validation of J Solutions for SC04 

The J solutions predicted by the SC04 model for internal and external surface cracks were 
validated against the results of manual calculations (see also Appendix 4). Excellent agreement 
was obtained between the two sets of computations. Comparisons were made when the applied 
load was specified in terms of an internal pressure and when the load was specified in terms of 
an arbitrary hoop stress distribution. In addition, internal consistency between the two forms of 
specifying the applied load was checked for two cylindrical geometries, corresponding to D/t 
equal to 22 and 102, where D is the outer diameter of the cylinder. In these cases, the arbitrary 
hoop stress was defined as that printed in the output when the load specification in terms of 
internal pressure is used. Figures 4 and 5 show the results obtained from this consistency check 
for internal and external cracks, respectively. It can be seen from the figures that there is 
excellent agreement between the pressure loaded solutions and the equivalent load defined in 
terms of an arbitrary hoop stress. 



Figure 4. Consistency between J estimations for internal cracks in cylinders (SC04) when the 
applied primary loading is specified in terms of an internal pressure and a hoop 
stress distribution corresponding to an internal pressure. 
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O a-tip (D/t=22) 
A c-tip (D/t-22) 
V a-tip (D/t=1 02) 
D c-tip (D/t=1 02) 
1 to t line 



J - SC04 - pressure loading 


Figure 5. Consistency between J estimations for external cracks in cylinders (SC04) when the 
applied primary loading is specified in terms of an internal pressure and a hoop 
stress distribution corresponding to an internal pressure. 


5.0 J SOLUTIONS FOR CORNER CRACKS IN PLATES SUBJECTED TO BENDING 
(Crack Model CC01) 

5.1 Implementation of J Solutions for CC01 (Bending) 

The current J solutions in NASGRO for comer cracks subjected to bending (CC01) are 
conservatively based. The accuracy of these solutions was improved by |performing elastic- 
plastic finite element analysis (FEA) to compute J solutions and use the results to reduce the 
conservatism in the J estimation technique used in the Phase 1 solutions in NASGRO. Only the 
main results of the FEA are presented in the main text of this Addendum, a more detailed 
description of the FEA modeling is provided in Appendix 1 . 

The results of the FEA were used to derive values of hi calculated as 


*,=• 


J-J. 




/ \ fl +1 

' M ' 


kKj 


(17) 


Values of hi were calculated for four different values of the strain-hardening exponent (n=l, 5, 
10, and 15) for each of the geometries in the analysis matrix. (The n= 1 values correspond to 
linear elastic solutions with a Poisson ratio equal to the plastic value of 0.5.) The resulting 
values of h] and the moment ratio ,(M/M‘ 0 ), at which they were evaluate are presented in Tables 
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2 through 5. The hi values were derived from the FEA J values at elliptical angles of 4.5° (c-tip) 
and 85.5° (a-tip) in order to avoid using the actual free surface values at 0° and 90° that are 
known to be subject to errors. The values of hi as a function of elliptical angle are shown 
graphically in Appendix 1. 

Table 2. Model CC01, Comer Crack, Bending Load, n=l. 


a/t 

a/c 

hi(a) at 85.5° 


0.2 

0.2 

0.7897 

0.3100 

0.2 

0.6 

0.5292 

0.5477 

0.2 

1.0 

0.3555 

0.5621 

0.5 

0.2 

1.1833 

0.8753 

0.5 

0.6 

0.4217 

1.2060 

0.5 

1.0 

0.2028 

1.1263 


Table 3. Model CC01, Comer Crack, Bending Load, n=5. 


a/t 

a/c 


hi(a)at85.5° 

hi(c)at45° 

0.2 1 

0.2 

2.145 

0.4408 

0.1719 

0.2 

0.6 

2.130 

0.3613 

0.2933 

0.2 

1.0 

2.122 

0.2396 

0.3367 

0.5 

0.2 

1.943 

0.7163 

0.4991 

0.5 

0.6 

2.316 

0.3014 

0.5051 

0.5 

1.0 

2.369 

0.1294 

0.4846 


Table 4. Model CC01, Comer Crack, Bending Load, n=10. 


ah 

a/c 


hi(a) at 85.5 ° 

hi (c) at 4.5° 

0.2 

0.2 

1.540 

0.2865 

0.1204 

0.2 

0.6 

1.508 

0.2568 

0.2125 

0.2 

1.0 

1.503 

0.1730 

0.2403 


0.2 

1.505 

0.3879 

0.3394 


0.6 

1.638 

0.1809 

0.2672 

0.5 

1.0 

1.652 

0.08035 

0.2488 


Table 5. Model CC01, Comer Crack, Bending Load, n=15. 


ah 

a/c 

KQB 

h 1 (a)at85.5° 

hj(c) at 4.5 ° 

0.2 


1.353 

0.19853 

0.08720 

0.2 

0.6 

1.343 

0.1843 

0.1551 

0.2 

1.0 

1.339 

0.1265 

0.1755 

0.5 

6.2 

1.380 

0.1945 

0.1971 

0.5 

0.6 

1.458 

0.09734 

0.1379 

0.5 

1.0 

1.463 

0.04455 

0.1279 
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The optimized yield moments and V values were determined for CC01 under bending using 
the FEA calculated hj values for n=l, 5, 10, and 15. As previously mentioned, the optimized 
scheme provides the values of V and M 0 that are independent of strain hardening exponent n and 

give the best fit between the RSM analytical approach and the FEA results. The results of 
applying this scheme are shown in Table 6 for the a-tip and c-tip as a function of a/t and a/c. 

Table 6. Optimized Yield Moments and V’s for CC01 under Bending. 


Comer Crack (CC01) Under Bending 



a-t 

ip 

c-tip 

a/t 

a/c 

Optimized 

4Mo*/t 2 Wcr 0 

V 

Optimized 

4Mo*/t 2 W<T 0 

V 

0.2 

0.2 

1.0068 

0.7224 

1.0055 

0.6310 

0.2 

0.6 

0.9931 

0.8745 

0.9896 

0.6666 

0.2 

1.0 

0.9877 

0.8576 

0.9879 

0.7399 

0.5 

0.2 

0.9607 

1.0072 

0.9237 

0.6641 

0.5 

0.6 

0.9436 

1.1279 

0.9618 

0.7002 

0.5 

1.0 

0.9379 

0.9840 

0.9658 

0.7313 



y Average 
Y a 

0.9289 

y Average 

0.6889 



(y Average + y Average yj = 0.8089 


In Figure 6, the optimized RSM results for hi (n>\) are shown plotted against the values derived 
from the FEA computations. The data points fall on or near the “1 to 1” line that represents 
100% accuracy for the optimized solutions, confirming that in principle the RSM approximation 
can attain high accuracy. As previously mentioned, the hybrid RSM solutions are employed in 
NASGRO in order to be able to determine J solutions for a wide range of a/t,\a/c , and h/c values. 
In Phase 1, these hybrid solutions do not use the optimized yield moments 'but instead use the 
expression for M* a given by equation (7). The values of hj predicted by the hybrid RSM 

solution for J are shown plotted against the FEA solutions in Figure 7. In this case, the hybrid 
RSM solutions consistently over-predict the FEA values, and the accuracy of the solutions is 
poor. This problem was attributed to the fact that equation (7) is not an accurate estimation for 
the optimized net section yield moment. With this in mind, studies were performed to obtain a 
modified form for M* 0 that increased the accuracy of the hybrid RSM solutions. The result of 

the investigation, the yield moment M ' 0 ccai , is shown in equation (18) where M * is given by 
equation (7). 


M' o CC0} =(1.033 + 0.184y)M o * (18) 

The NASGRO EPFM CC01 bend solution for J p in Phase 3 is taken, therefore, as the hybrid 
RSM solution given by the equation 
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*CC01 
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V* 



(19) 


where the value of V is 0.8089 (see Table 5), and The corresponding hybrid RSM solutions 
for hi(n) are given by 


h l (n) = 


J t HV 


<wi 


M 


K M o 


♦CC01 


( 20 ) 


The predictions of equation (20) for n>l are compared to the FEA solutions for the a-tip and 
c-tip obtained from Tables 3 through 5 in Figure 8. It can be seen that there is a significant 
increase in the accuracy of the hybrid RSM solutions obtained using M* CC01 compared to those 

solutions using Af* . Indeed, the new solutions are evenly scattered about the 1 to 1 line rather 
than consistently over-estimating the values of hj(n). 



Figure 6. Comparison of hj(n>l) for CC01 (bending) computed using FEA with the results 
obtained from applying the RSM using optimized net section yield moments. 
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Figure 7. Comparison of h](n>l) for CC01 (bending) computed using FEA with the results 
obtained from applying the hybrid RSM using equation (7) for the net section yield 
moment. 



Figure 8. Comparison of hj(n>l) for CC01 (bending) computed using FEA with the results 
obtained from applying the hybrid RSM using equation (18) for the net section 
yield moment. 


16 





5.2 Validation of J Solutions for CC01 (Bending) 

The CC01 J solutions were implemented for bending using the existing CC01 SDF routines, 
the net section yield moment given by equation (18) and the average V value of 0.8089 given in 
Table 6. 

Figure 8 provides validation for the new solutions against FEA results in terms of the hybrid 
RSM and FEA solutions for the hi functions, defined according to equation (20). This figure 
demonstrates the kind of accuracy that can be obtained from the RSM solutions. Additional 
validation was obtained by comparing manual calculations for J p with J p values computed using 
the EPFM module. These computations confirmed that the hybrid RSM solution for J had been 
accurately implemented in the computer code (see Appendix 4). 

6.0 J SOLUTIONS FOR EMBEDDED CRACKS IN PLATES SUBJECTED TO 
TENSION AND ARBITRARY STRESSES (Crack Models EC01/EC02) 

6.1 Implementation of J Solutions for EC01/EC02 

The Phase 1 J solutions in NASGRO for embedded cracks subjected to tension (EC01) were 
considered conservatively based. (However, as will be seen below, this proved not to be the 
case.) As was done for CC01 in bending, the accuracy of these Phase 1 solutions was improved 
using the results of FEA to compute J solutions. Only the main results of the FEA are presented 
in the main text of this Addendum, a more detailed description of the FEA modeling is provided 
in Appendix 1. 

The results of the FEA were used to derive values of hi calculated as 


*»=■ 


J-J. 


™o £ o 


\ n+1 


kKj 


( 21 ) 


Values of hj were calculated for four different values of the strain-hardening exponent (n=l, 5, 
10, and 15) and the resulting values of hi and the load ratios (P/P* ) at which they were 
determined are presented in Tables 7 through 10. The hi values were derived from the FEA J 
values at elliptical angles of 4.5° (c-tip) and 85.5° (a-tip) in order to avoid possible errors at 0° 
and 90°. The values of h] as a function of elliptical angle are shown graphically in Appendix 1. 

Table 7. Embedded Crack, Tension Load, n=l. 


2aJt 

a/c 

hi(a) at 85.5 ° 

hi(c) at 4.5 ° 

0.2 

0.2 

0.4337 

0.09056 

0.2 

0.6 

0.2930 

0.1755 

0.2 

1.0 

0.1959 

0.1959 

0.5 

0.2 

1.3037 

0.2442 

0.5 

0.6 

0.8013 

0.4685 

0.5 

1.0 

0.5119 

0.5062 


17 























Table 8. Embedded Crack, Tension Load, n=5. 


2at t 

a/c 


hi(a) at 85.5 ° 

hj(c) at 4.5° 

0.2 

0.2 

2.134 

0.8231 

0.1851 

0.2 

0.6 

2.129 

0.5299 

0.3247 

0.2 

1.0 

2.212 

0.3662 

0.3654 

0.5 

0.2 

2.354 

4.3138 

0.7913 

0.5 

0.6 

2.403 

1.9887 

1.1550 

0.5 

1.0 

2.329 

1.1850 

1.1319 


Table 9. Embedded Crack, Tension Load, n=10. 


2aft 

a/c 

WEBm 

hi(a) at 85.5° 

hi (c) at 4.5° 

0.2 

0.2 

1.521 

1.0699 

0.2376 

0.2 

0.6 

1.458 

0.6843 

0.4087 

0.2 

1.0 

1.515 

0.4755 

0.4732 

0.5 

0.2 

1.606 

8.8905 

1.7439 

0.5 

0.6 

1.619 

3.4996 

2.0009 

0.5 

1.0 

1.589 

1.9506 

1.8448 


Table 10. Embedded Crack, Tension Load, n=15. 


2a/t 

a/c 

BM 

hj(a) at 85.5° 

hi (c) at 4.5° 

0.2 

0.2 

1.340 

1.2678 

1 0.2808 

0.2 

0.6 

1.284 

0.8071 

i 0.4765 

0.2 

1.0 

1.333 

0.5613 

0.5569 

0.5 

0.2 

1.409 

14.800 

3.1478 

0.5 

0.6 

1.413 

5.3330 

3.0537 

0.5 

1.0 

1.396 

2.8317 

2.6808 


The optimized yield loads P* and V values were determined for EC01/EC02 using the FEA 

calculated hi values for n= 1, 5, 10, and 15. The results of applying this scheme are shown in 
Table 11 for the a-tip and c-tip as a function of aJt and a/c. 
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Table 11. Optimized Yield Moments and V’s for EC01/EC02 under Tension. 


Embedded Crack (EC01/EC02) Under Tension I 



a-tip 

c-tip 

2aJt 

a/c 

Optimized 

Po*/tW<T 0 

V 

Optimized 
Po*/tWi T 0 

V 

0.2 

0.2 

0.9181 

1.5916 

0.9194 

1.7213 

0.2 

0.6 

0.9191 

1.5229 

0.9226 

1.5788 

0.2 

1.0 

0.9181 

1.5656 

0.9186 

1.5660 

0.5 

0.2 

0.7925 

1.9930 

0.7812 

1.8418 

0.5 

0.6 

HEXIEiH 

1.6569 

0.8156 

1.6609 

0.5 

1.0 


1.6202 

0.8246 

1.5714 



y Average 
v a 

1.6584 

y Average 
V c 

1.6567 



( yMcragr + yAverag, )/2 = j 5575 


In Figure 9, the optimized RSM results for hj(n>l) are shown plotted against the values derived 
from the FEA computations for EC01. The data points fall on or near the “1 to 1” line that 
represents 100% accuracy for the optimized solutions, yet again confirming that the RSM 
approximation can attain high accuracy. 



hjn) derived from FEA 


Figure 9. Comparison of hi(n>l ) for EC01/EC02 computed using FEA with the results 

obtained from applying the RSM using optimized net section yield loads. 


The hybrid RSM J solutions for EC01 in Phase 1 are generated using the equation 
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(22) 


jRSM 

J P 


= J e (a,c)juVa\ 


( p V 1 


\ P 0 J 


where P* is related to the reduced load bearing area of the plate and is given by 

P;=<r 0 (tW-mc) (23) 

The values of hi corresponding to this hybrid RSM solution for J are shown plotted against 
the FEA solutions in Figure 10. In this case, the hybrid RSM solutions consistently under- 
predict the FEA values, and the accuracy of the solutions is poor. This problem is due to the fact 
that equation (23) provides a poor representation of the optimized net section yield moment. 



Figure 10. Comparison of hj(n>l ) for EC01/EC02 computed using FEA with the results 
obtained from applying the RSM using equation (23) as the net section yield load. 

Studies were performed to obtain a modified form for P* that increased the accuracy of the 

hybrid RSM solutions. Based on this investigation, the Phase 3 NASGRO EPFM EC01 and 
EC02 tension solutions for J p are taken as the hybrid RSM solution: 


f 


v -i 


r RSM 
J P 


_ ,2a a 2c. tr , 

7 - ( TV¥^ 


p*EC01/ ECOlf 2fl 2c 

0 u Vw, 


(24) 
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1-v* 

In this equation, the value of V is taken as 1.6575 (see Table 11) and p = -r for both the 

1-v; 


a-tip and c-tip, where subscripts e and p signify the elastic and plastic values, respectively, of 
Poisson’s ratio. The net section yield load is given by: 


p'ECOMECOi _ /I A 4 2d 2c .| 
» V t W n 


1 + 0.025 


V 


1 + 401-1 

W W. 


\U/ \^/2^ 




(, _2a2c"\ 

, 2a 2c (a V t l T7 f] J 

1-0.9 

1 + 0.3102 v ' 

t wj 

t W\c) 


\ y 


{ 2a 2c Y (a 

t w) U 


P' alc< 1 


(25a) 


P* £C0,,£C02 =(1-0.4— — )| 
* t JV 1 


f 


2a2c('cY s ( , -l < TFn 


V 


fl-0.9— — 
t w , 


1 + 0.3102 


t W 


a) 


P* a/c> 1 


(25b) 


where P* is given by equation (23). The hybrid RSM solutions for hj(n) are given by 


h 1 (n) = - 


J e (a,c) M V 


a ° S °2 


V 


(26) 


r>*£C01/£C02 

V - * 0 J 


The predictions of equation (26) for «>1 are compared to the FEA solutions from Tables 8 
through 10 for the a-tip and c-tip in Figure 11. The accuracy of the solutions is greatly improved 
using the modified form for the net section yield load, and the results are now evenly scattered 
about the 1 to 1 line. 
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h^n) derived from FEA 

Figure 11. Comparison of hj(n>l) for EC01/EC02 computed using FEA with the results 
obtained from applying the RSM using equation (25) as the net section yield load. 

In the EC02 model the applied load is not explicitly defined. Instead, the user specifies an 
arbitrary stress distribution, and the tensile force, P, used in equation (**) is obtained by 
integrating this stress distribution. 

6.2 Validation of J Solutions for EC01/EC02 

The new FEA J solutions for embedded cracks subjected to tension loading were used to 
update the V values and net section yield solutions in the Phase 1 EPpM module, and to 
implement the new EC02 solutions generated in Phase 2. The SIF solutions' employed in EC02 
are based on the solutions in KCALC, a program for computing SBFs for cracks in arbitrary 
stress fields developed and copyrighted by Southwest Research Institute® (SwRI®). This 
program was used because, unlike the case for SC02, NASGRO did not have the capability of 
calculating SDFs for embedded cracks subject to arbitrary stress fields. KCALC routines have 
been validated and are employed in several programs developed by SwRI, such as DARWIN™ 
(Design Assessment of Reliability With INspection), a software design code, developed for the 
Federal Aviation Administration (FAA) to help engine manufacturers improve the safety of jet 
engines used in commercial airliners. 

The new J solutions were partly validated by comparing the EC01 and EC02 solutions for 
uniform stressing. The results are shown in Figure 12 where J values derived from the EC01 
model are plotted against J values computed using the new EC02 model solutions. Perfect 
agreement between the two sets of solutions occurs when the data points fall on the 1 to 1 line. It 
can be seen from Figure 12 that excellent agreement obtains between the EC02 and EC01 
solutions, indicating that the integration routine used to determine the external force from the 
arbitrary stress distribution specified in EC02 and the resulting net section yield solution is 
correctly calculated. The small differences between the EC01 and EC02 J values arise because 
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the NASGRO SIF solution for uniform stressing is used in EC01, whereas, as previously 
mentioned, the KCALC SIF solution is used in EC02. 



j • EC01 


Figure 12. Comparison of J estimations obtained using EC01 and EC02. The arbitrary stress 
fields used in the EC02 computations were chosen to simulate uniform tension. 


An additional verification test for the EC02 model geometry was performed. This was based 
on an independently developed computer program that employed the KCALC routine and the 
same net section yield load equations as used in EC02. This independent program was used to 
generate J values (hereafter referred to as J estimation values) against which the NASGRO 
EPFM EC02 solutions could be compared. In this comparison, two different load cases were 
used. The first consisted of a primary load corresponding to a linear stress field of the form 

140 - 80— ksi. The second load case involved combined primary and secondary loads, with the 


stress distribution for the primary load given by a uniform stress equal to 100 ksi, and the 
secondary load corresponding to a self-equilibrated stress of the form 

,2 


-200 + 1200 


T'Kt) 


ksi. The J values determined using NASGRO for these two load 


cases are plotted against the J estimation values in Figures 13 and 14. Agreement between the 
two sets of solutions occurs when the data points fall on the 1 to 1 lines shown in the figures. It 
can be seen that excellent agreement is obtained between the NASGRO routine solutions and 
those obtained using the independently developed program. 


Additional validation for EC02 J solutions based on manual calculations is provided in 
Appendix 4. 
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J - EC02 estimation 

Figure 13. Comparison of EPFM Module J solutions for EC02 with independently derived 
solutions (J estimates) that used KCALC SIF solutions, and the same V and net 
section yield loads used in the Module. The primary load is represented by a linear 
stress distribution. 



Figure 14. Comparison of EPFM Module J solutions for EC02 with independently derived 
solutions (J estimates) that used KCALC SIF solutions, and the same V and net 
section yield loads used in the Module. The primary load corresponds to a uniform 
stress and the secondary load is a self-equilibrated quadratic stress. 
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7.0 TECHNICAL ISSUES RELATED TO PHASE 3 ENHANCEMENTS 


The major technical issues needed to be overcome to implement the proof test analysis 
modules were the development of: 2-DOF failure analysis routines for critical load and critical 
crack analyses; tear-fatigue routines that accurately included 2-DOF interactions between the 
a-tip and c-tip; and reliability analyses for MCPT analysis. 

7.1 Ductile Failure Analysis Routines for 2-DOF Cracks 

The development of 2-DOF ductile failure analysis routines for critical crack and critical load 
analyses is a major advance on the 1-DOF failure routines incorporated into Phase 1 and 2, and is 
a necessary enhancement in preparation for the introduction of proof test analysis modules in 
NASGRO. 


The conditions for ductile instability at the a-tip for 1-DOF cracks are defined by the 
equations: 


J a (®| » P) - Jr (Afl, ) 


(27) 


dJ m _ dj* 

da d(Aa t ) 


where a, is the initial crack depth before tearing occurs and Aa, is the amount of tearing at 
instability. These equations state that instability will occur when the applied J equals Jr and 
simultaneously the J curve is tangential to the Jr curve. 

For 2-DOF cracks, these conditions become: 

J a (a i + A a , , + Ac, ,P) = J R (A a , ) 


J c (ai +Aa,,c ( . + Ac t ,P) = J R (Ac t ) 


(28) 


. dc 


a,P V 


K 

da 


I 

K.p <*(*“.) 


(ay w dj, 

• 3a } < F <i(Ac,) 



where a, and c, are the initial crack depth and initial half surface length, respectively, and Ac, is 
the amount of tearing at the c-tip at instability. 

The 2-DOF instability conditions show that instability does not occur when a 1-DOF 
instability condition occurs at either the a-tip or the c-tip, but that instability is dependent on the 
conditions at both the a-tip and c-tip and does not correspond to a tangency point, as does the 
1-DOF case. Indeed, the instability condition for cracks with 2-DOF states that instability will 
only occur when both the a-tip and the c-tip are simultaneously unstable. 
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7.2 Proof Test Module 


The proof test module is based on the NASA Final Report, “ Guidelines for Proof Test 
Analysis, ” delivered to MSFC under NASA Contract NAS8-39380. Reference should be made 
to this document for more details concerning proof test design and analysis. Herein, only a brief 
summary of the proof test modules developed in Phase 3 is given. 

Figure 15 provides an overview of the routines included in the proof test module. Two types 
of analysis can be performed as part of the proof test procedures: either a Safe Life Analysis or a 
Proof Test Analysis. Two options are available if the Safe Life Analysis is selected: either 
Critical Flaw Size or Fatigue Life. The purpose of the Safe Life Analysis is to enable the proof 
test analyst to perform a pre-proof test calculation to determine those regions of a component 
that may be life limiting. The proof test should be designed to screen out unacceptable flaws in 
these regions. The life limiting regions may be defined in terms of low fracture tolerance for 
small cracks or in terms of low fatigue life. High stresses and/or low toughness may give rise to 
low flaw tolerance, and high cyclic stress ranges and/or environmental factors may give rise to 
fast crack propagation rates and low fatigue lives. 

Three options are currently available if the proof test analysis option is selected: either Proof 
Load Analysis or Flaw Screening Analysis or Final Crack Sizes. The purpose of the proof load 
analysis is to determine the proof load necessary to screen out flaws above a specified size. The 
purpose of the flaw screening analysis is to determine the flaw sizes that are screened out by a 
specified proof load. The final crack sizes option enables analysts to determine the increase in 
sizes of specified flaws due to application of the proof load. Although implementation of a proof 
load analysis or a flaw screening analysis will provide an analyst with information regarding 
which sizes of flaws will not be present in the component after it has been proof tested, these 
options will not predict how the population of flaws that survive the proof test has grown due to 
ductile tearing that did not result in crack instability. The final crack size dption is intended to 
provide this information. ' 
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Proof Test Procedures 
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Figure 15. Overview of the Proof Test Module incorporated into NASGRO. 





7.3 Tear-Fatigue 

Tear-fatigue occurs when load cycling is severe enough to result in simultaneous fatigue 
crack extension and ductile tearing. The synergy between these two mechanisms of crack 
propagation results in enhanced crack propagation rates with respect to fatigue crack growth. 
This point is illustrated in Figure 16 which shows measured crack growth rates plotted against 
the applied closure corrected cyclic SIF, AK e ff, for stress ratios, R, of 0.5, 0.1, and -1. These 
results demonstrate that the enhanced crack growth rate due to tear-fatigue can be an order of 
magnitude higher than the predicted fatigue crack growth rate, that tear-fatigue can occur at any 
R value and, unlike fatigue crack growth, that AK e jj does not collapse the growth rate in the tear- 
fatigue regime onto a single curve. A pictorial representation of the tear-fatigue process is 
shown in Figure 17 and illustrates how the mechanism depends on both fatigue crack growth 
properties and the J-R curve of the material characterizing the resistance to tearing. 

The tear-fatigue methodology is only applicable to ductile materials and is needed to 
implement the MCPT module in NASGRO. The tear-fatigue methodology is limited to cases of 
single amplitude loading, which makes it suitable for applying to an MCPT analysis where the 
proof load is repeatedly applied and removed. 

There are two stages to implementing tear-fatigue for 2-DOF flaws. The first stage consists 
of calculating the amount of ductile tearing that occurs on first application of the proof load. 
Since the crack tip driving forces at the a-tip and c-tip change as tearing occurs, this stage 
involves incrementally increasing the applied load up to its maximum value in the fatigue cycle 
taking into account the resulting incremental changes in the tear lengths at the two tip positions. 
The second stage consists of actual tear-fatigue as the applied load is cyclically applied, and the 
a-tip and c-tip incrementally increase in length after each cycle due to fatigue crack growth and 
ductile tearing. 

The two stages can be expressed mathematically the following equations, j 
Stage 1: First (monotonies load application 


The incremental changes in the tear lengths at the two tips, 5a t and 5c t due to an incremental 
change in applied load, SP are given by: 


Sa t = SP — - 

(j - r c ' p> i( t - j a ' p i c ' p 

W RM, J a,a)\ J RM l J c,c ) J a,c J c,a 


(29) 


j a < c ( j _ J c ’ p \+ J e ' r T a ' c 
^ c,P\ J RM, J a,a J c,a J a,P 


t c t P ra t c 


5c t = SP- 


Vrm, -jHWrm, -W-kp:: 




ja,P jc,P 
:,a 


(30) 


In these equations, subscript t refers to tear, a and c refer to the a-tip and c-tip, respectively, P to 
the applied load. Also, the following abbreviations are used: 
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( 31 ) 


J*L] 

x,y UJ, 


RM, 


_Jh L 


dJ. 


4 Aa <) 


R,Ac, 


d(Ac t ) 


(32) 


where A a, and Ac, are the current tear lengths at the a-tip and c-tip, respectively, and Jr is the 
resistance curve. 


Stage 2: Tear- fatigue (cyclic loading} 


The incremental changes in crack lengths at the two tips due to a single fatigue cycle are 
given by: 

da f = A{U a<eff Y 


dc f = 


(33) 


In these equations, subscript / refers to fatigue, subscript effio a crack closure corrected quantity, 
and A J a tff and tj a tf[ are the cyclic changes in /at the a-tip and c-tip, respectively. 

The corresponding incremental changes in the tear lengths are assumed to occur at maximum 
load, Pmax, in the cycle and are given by: 


a _ +■#-&/) 


\ J RM, J a,a )\ J R,6c, J c,c ) J a,e J c. 


‘fm 


r“’ P m 


jC,P m 

:,a 


(34) 


& _ v£-sc f 


(35) 


Equations (29), (30), (34), and (35) show that conditions for instability occur when the 
denominators in these equations become zero, and that these conditions are the same for 
monotonic loading and cyclic tear-fatigue and for both the a-tip and the c-tip. These instability 
conditions are precisely those specified in equation (28). 
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Figure 16. Measured crack growth rate data showing how tear-fatigue accelerates the growth 
rate with respect to fatigue, illustrated by the Paris equation fit to the data. 



• before fatigue cycle the 
crack is stable and J=J R 

• crack advances by fatigue 

• resistance curve process 
zone moves with crack tip 

• applied J increases and 
exceeds resistance (J>J R ) 

• crack tip tears until stable 

(J=Jr) 

• crack advances by fatigue 
on next cycle etc 


Figure 17. Illustration of how the mechanism of tear-fatigue involves synergy between fatigue 
crack growth and ductile tearing. 
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7.4 Multiple Cycle Proof Test Analysis (MCPT) 

Normally, a component is subjected to a single load cycle during proof testing before 
entering service. However, for ductile materials, it has been observed that applying multiple 
load cycles can increase service reliability in some circumstances compared to a single cycle 
proof test. 

According to deterministic proof test analyses, MCPT will cause flaws to extend so a 
component will enter service with a larger flaw size population than would be the case without 
MCPT, reducing service reliability. A probabilistic calculation is needed to demonstrate that 
MCPT can increase service reliability. The argument is based on the fact that MCPT will 
beneficially change the service reliability by removing those components with large flaws that 
are service life limiting before they enter service, more than compensating for the potential 
increase in flaw size population in those components that survive the MCPT. 

The methodology employed in the NASGRO MCPT module is based on the work described 
in the Final Report "A Comparison of Single Cycle Versus Multiple Cycle Proof Test Strategies ” 
performed under Contract NAS8-37451. Consistent with that methodology, there is only one 
random variable considered in the probabilistic analysis, namely, the initial crack depth. 

The MCPT module calculates the following failure probabilities: probability of failure for 
proof cycles only; probability of failure for proof plus service cycles; and the conditional 
probability of failure in Ns' service cycles given no failure in Np 'proof cycles. 

The reduction of the problem to a single random variable (crack size) allows the probability 
problem to be reformulated in terms of initial crack size. Therefore, the probability of the 
number of service cycles being less than or equal to a prescribed number of service cycles is 
expressed mathematically as 

?[^^Afe'] = P[//(a f )^A&'] = P[fl < ^/r , (^j , )] = P[a^a;] (36) 

where Ns is the number of service cycles at failure, Ns 'is a specified number of service cycles, a, 
is the initial crack size random variable. Ns = H(aj) denotes the crack growth function, H~‘ is the 
inverse of the crack growth function, and a\ is the initial crack size that causes failure on the Ns ' 

service cycle. Similarly, the probability of the number of proof cycles (N p ) being greater than a 
prescribed number of proof cycles (Np ) is expressed mathematically as 

P[Np > Wp'] = />[/%)> NP'] = P[°< < H-' (Np')] = P\a, < a'] (37) 

where Np is the number of proof cycles at failure, Np 'is a prescribed number of proof cycles, af 

is the initial crack size which causes failure on the Np '+1 proof cycle. The initial crack sizes for 
both Np ' proof cycles and Np ' (proof) + Ns ' (service) cycles are printed in the output file along 
with the probabilities defined in equations 1 and 2. The final probability value calculated is the 
conditional probability of failure in Ns ' service cycles given no failure in Np ' proof cycles. 
Mathematically, the conditional probability is expressed as 


31 



r , , , -i P [ (Ns < Ns) n (Np>Np) 

p [(***) |w>j»)]- i J 


( 38 ) 


Because both P [ Ns < Ns'] and P [ Np > Np] can be represented in terms of the initial crack 
size distribution a h the intersection term in equation (3) can be computed algebraically, once the 
initial crack sizes af and a ‘ are known. 

An outline of the approach is shown in Figure 1 8 that illustrates the three stages involved. 


The MCPT module currently only 
has one random variable, the initial 
flaw size before the proof test. This 
enables analytical simplifications to 
be made, avoiding lengthy Monte 
Carlo-calculations. 


Stage 3: The probability of the 
component failing within the 
specified service lifetime given that 
it survives the MCPT is the cross 
hatched area shown in the figure. 


Stagel : Determine initial flaw 
size to just survive the MCPT plus 
the specified service life 




probability distribution function 
for initial flaw sizes 


Stage 2: Determine the 
initial flaw size to just 
survive the MCPT 


Figure 18. The three stages in the conditional probability calculations in the MCPT module. 


The MCPT module can be applied to either 1 -degree of freedom (DOF) flaws, or 2-DOF 
flaws. In the case of the latter, tear-fatigue crack growth is calculated using 2-DOF, based on the 
monotonic and cyclic crack tip driving forces at the deepest and surface points. The service 
lifetime calculation begins at the end of the MCPT and uses the final crack size at the end of the 
proof test as the initial size. The service lifetime calculations are again based on 2-DOF crack 
growth routines if the problem involves 2-DOF flaws. Note that tear-fatigue crack growth is not 
allowed for under service conditions as, in general, these will involve variable amplitude loading 
for which the tear-fatigue routines are not applicable. 

The MCPT module calculates the conditional probability of failure for a component for a 
user specified service lifetime given that the component survives the MCPT. The MCPT is 
advantageous if this probability is less than the probability of failure determined for a single 
proof test cycle, or when no proof test is performed. 
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Although the probabilistic part of the calculations are performed analytically rather than 
employing Monte Carlo or other numerical methods, never-the-less significant computation time 
is needed to search and find the initial crack depths that will grow to failure in the user specified 
proof and service cycles, especially for 2-DOF flaws. Thus, the calculations to evaluate failure 
probability for a single pair of user specified values for the number of proof cycles and service 
cycles may take several minutes or more, depending on the speed of the computer used. 

8.0 EXAMPLES, VALIDATION, AND PROGRAM ISSUES 

8.1 Example Input and Output for Running the EPFM/Proof Test Modules 

Examples of the input data needed to run the Modules are presented in Appendix 2. This 
Appendix contains ten tables listing the data necessary to interactively input data to create the ten 
example files, Examl.inp through ExamlO.inp, contained on the distribution CD. The data is 
presented in the order requested by the screen prompts from the Modules. The tables list the 
name of the input parameter, its value, the units of the parameter, and a brief description of it. 
The ten examples are summarized in Table 12. 


Table 12. Summary of analyses performed in Examples 1 through 10. 


Example 

No. 

Option 

No. 

Analysis Type 

Crack Model 

Loading 

1 

5 

J computation 

SC04 


2 

6 

critical crack 

SC02 


3 

6 

critical load 

EC02 

primary 

4 

7 

fatigue life 

CC01 


5 

8 

pre-proof test critical 
crack 

SC02 

primary and secondary 

6 

8 

proof test: flaw 
screening 

SC04 

primary 

7 

8 

proof test: proof load 

EC02 


8 

8 

proof test: final crack 
size 

CC01 

primary 

9 

9 

tear-fatigue 

SC04 


10 

10 

MCPT 

SC02 

primary 


The output files, Examl.out through ExamlO.out, respectively, produced by the example 
input files are also contained on the distribution CD. 

Hard copies of the ten input files and the corresponding output files are given in Appendix 3. 

8.2 Validation 

The validation of the Modules has been largely directed at ductile failure and fatigue 
analyses, and, in particular, those analyses that involve tearing and tear-fatigue with 2-DOF. As 
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mentioned previously, the 2-DOF calculations for ductile materials proved the most difficult to 
computationally implement. 

Appendix 4 lists in tabular form the results of part of the exercise performed to validate the 
Modules. This exercise complements and provides additional verification material to that 
already presented in the main part of this Addendum. Appendix 4 presents validation for all the 
EPFM options (Options 5 through 10). Except for the Option 5 (J calculations), all the analyses 
used in the validation involved 2-DOF, and except for Option 5 and Option 7 (fatigue lifetime), 
all the validation analyses addressed ductile fracture behavior. 

Appendix 4 presents the results of applying two methods for validating the Modules. In the 
first, manual spreadsheet calculations were performed to independently evaluate the results of 
applying the Modules. These validation runs are summarized in Tables A4.1 through A4.6, and 
Table A4.15 in Appendix 4. The verification runs performed in this exercise are listed in 
Table 13. 


Table 13. List of the manual spreadsheet calculations performed to validate the NASGRO 
Modules. More details are given in Appendix 4 to which the table numbers refer. 


Table 

Number 

Crack 

Model 

Option 

Number 

Description 

A4.1 

SC04 

5 

J estimation 

Module results for Example 1 in Appendix 2 verified against 
manual spreadsheet calculations by comparing predicted J p 
values. 

A4.2 

SC02 

5 

J estimation 

Module results verified against manual spreadsheet calculations 
by comparing predicted J p values. Primary and secondary 
loads. Primary stress distribution integrates to a tensile force 
and zero moment. 1 

A4.3 

SC02 

5 

J estimation 

Module results verified against manual spreadsheet calculations 
by comparing predicted J p values. Primary stress distribution 
integrates to a tensile force and moment. 

A4.4 

CC01 

5 

J estimation 

Module results verified against manual spreadsheet calculations 
by comparing predicted J„ values. Primary bending load. 

A4.5 

EC02 

5 

J estimation 

Module results verified against manual spreadsheet calculations 
by comparing predicted J p values. Primary and secondary 
loads. 

A4.6 

SC02 

6 

Critical crack 
size 

Module results for Example 2 verified against manual 
spreadsheet calculations based on running Option 5 to obtain J 
estimates. Results demonstrate that that the applied J values at 
the a-tip and c-tip fall on the J-R resistance curve, and the 
ductile instability criterion is satisfied. 

A4.15 

SC02 

10 

MCPT 

Module results for Example 10 verified against manual 
spreadsheet calculations by comparing predicted conditional 
probability of failure value. The probabilities are evaluated 
using the initial crack sizes calculated by the module for cracks 
that would just survive the proof test and service lifetime, 
respectively. 
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In the second verification method, self-consistency checks were performed for the Modules 
by calculating the same result twice using different options and showing that similar results were 
produced. The results of these internal consistency checks are summarized in Tables A4.7 
through A4.14, and Tables A4.16 and A4.17 in Appendix 4. The verification runs performed are 
listed in Table 14. 

Table 14. List of the internal consistency calculations performed to validate the NASGRO 
Modules. Self-consistency between the Modules is investigated by using two 
different options to calculate the results for similar problems. More details are 
given in Appendix 4 to which the table numbers refer. 


Table 

Number 


Crack 

Model 


Options Compared 


Option 


Option 


Description 


SC02 


6 

Critical 

Crack 


6 

Critical 

Load 


The critical crack sizes (Option 6) determined in 
Example 2 are used to specify the initial crack 
sizes in critical load (Option 6) calculations. 
Internal consistency is achieved by demon- 
strating that the critical load equals the applied 
load used in the critical crack size computations, 
and the predicted tear lengths are the same. 


EC02 


6 

Critical 

Load 


6 

Critical 

Crack 


The critical load (Option 6) results calculated in 
Example 3 are used to specify the applied loads 
in critical crack size (Option 6) calculations. 
Internal consistency is achieved by demon- 
strating that the critical crack sizes equal the 
initial crack sizes used in the critical load 
computations, and the predicted tear lengths are 
the same. 


CC01 


7 

Fatigue 

Life 


8 


Safe Life: 
Fatigue 
Life 


Fatigue crack growth behavior predicted in 
Example 4 using the fatigue life analysis (Option 
7) is shown to be consistent with similar 
behavior predicted by the Safe Life: Fatigue Life 
analysis (Option 8), 


10 


SC02 


8 

Safe Life: 
Critical 
Crack 
Size 


6 

Critical 

Load 


The critical crack sizes predicted by Example 5 
in a Safe Life analysis (Option 8) are used to 
specify the initial crack sizes in critical load 
(Option 6) calculations. Internal consistency is 
achieved by demonstrating that the critical load 
equals the applied load used in the critical crack 
size computations, and the predicted tear lengths 
are the same. 


11 


SC04 


8 

Proof 
Test: 
Flaw 
Screening 


8 

Proof 

Test: 

Proof 

Load 


The screened crack sizes predicted in a Proof 
Test analysis (Option 8) in Example 6 are used to 
specify the initial crack sizes in a proof load 
estimate using Proof Test analysis (Option 8). 


Internal consistency is achieved by 











Table Crack I Options Compared 
Number Model 



Description 


demonstrating that the proof load needed to 
screen against the initial crack sizes equals the 
proof load used in the flaw screening 
computations, and the predicted tear lengths are 
the same. 

The screened crack sizes predicted in a Proof 
Test analysis (Option 8) are used to specify the 
initial crack sizes in a critical load analysis 
(Option 6). Internal consistency is achieved by 
demonstrating that the critical load corresponds 
to the proof load used for flaw screening. 

The screened flaw sizes proof load used in the 
Proof Test analysis (Option 8) are demonstrated 
equal to the instability crack sizes calculated 
using a critical crack analysis (Option 6). used to 
specify the initial crack sizes in a critical load 
analysis (Option 6). 


In Example 7 the proof load (Option 8) 
necessary to screen against a specified initial 
crack size is shown to be consistent with the 
predicted flaw size screened against when this 
load is applied as the proof load. 


The final crack size at the end of a Proof Test 
(Option 8) in Example 8 is demonstrated to be 
the same as the final crack size at the end of the 
first load in a Tear-Fatigue (Option 9) analysis. 


The final crack size at the end of the first load 
application in the Tear-Fatigue (Option 9) 
analysis of Example 9 is demonstrated to be the 
same as the final crack size at the end of a Proof 
Test (Option 8). 


The calculated initial crack size for a specified 
service life in a MCPT (Option 10) analysis with 
no proof test is shown to be consistent with a 
fatigue life (Option 7) analysis. 


The calculated initial crack size for a specified 
number of proof test cycles in a MCPT (Option 
10) analysis with no service cycles is shown to 
be consistent with a Tear-Fatigue (Option 9) 
analysis. 
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8.3 Program Issues 

8.3.1 Problems Occurring During Program Execution 

In some instances the search routines employed to solve Option 6, Option 8, and Option 10 
analyses may encounter problems. The causes of these problems are usually attributable to one 
of the following: 

1. Limitations on the range of geometrical parameters (a/c, a/t) for which the Phase 1 
SIF solutions are valid. If the required solution falls outside of these validity ranges 
then the search routines in the program will fail. 

2. Critical crack size and critical load calculations for ductile materials involving 2-DOF 
cracks involve evaluating derivatives of J and Jr (for example, see equation (28)). 
The search routines may encounter problems in finding solutions in these cases 
because of discontinuities in the derivatives caused by: 

A. The change in gradient in Jr as the J-R curve transitions from the blunting line 
to the ductile tearing curve; 

B. The change in gradient in Jr as the J-R curve transitions from the ductile 
tearing curve to the saturation value where the gradient becomes zero; 

C. The transition of J from a continuously varying function of crack size and 
load to an assumed infinite value when the reference stress equals or exceeds 
the flow stress defined in equation (3). 

However, the user is recommended to check the reasonableness of input data before first 
assuming that program problems are caused by one of the reasons given above. 

8.3.2 CPU Time 

In some cases, the number of iterations needed to accurately compute the 2-DOF results for 
ductile materials is very large due to the sensitivity of the results to growth history. The 
computations are particularly long if both primary and secondary loads are applied. As a result, 
in deterministic calculations the computations can take between seconds to tens of seconds to 
complete. The root finding procedures necessary to implement the MCPT analysis involve even 
longer computations and, in these cases, CPU times that may extend out to minutes in duration. 

8.3.3 CD Contents 

The delivered CD contains the following items: 

1 . An electronic version of this Letter Report. 

2. An executable file for running the MSFC Version 6.0 of the NASGRO EPFM and Proof 
Test Modules. 

3. Input files for exercising the executable and the corresponding output files. 
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APPENDIX 1: Finite Element Analysis of CC01 (Bending) and EC01/EC02 (Tension) 


A.l: FEA for CC01 (Bending) Model 

Finite element models were created for CC01 geometries. A schematic of the CC01 
model is shown in Figure A.l. Each of the CC01 finite element meshes was generated for the 
present matrix of crack geometries using Patran. The finite element modeling took advantage of 
appropriate symmetry conditions to reduce the size of the models needed for analysis. Thus, in 
the case of CC01, symmetry conditions enabled the model size to be reduced to half the size 
needed to model the full geometry. Consistent with the FEA-based J results used in Appendix K, 
the ratios b/c=4 and c/h=0.25 were held constant for all the analyses, where h is the height of the 
cracked plate. The elements used in the analysis were 20-noded brick elements with reduced 
integration. The 20-noded brick elements utilized quadratic shape functions for improved 
accuracy under bending conditions. Additionally, the reduced integration element enabled more 
accurate representation of the constant volume condition associated with plastic deformation. 
Each finite element model contained a focused ring of element around the crack front. Crack tip 

n_ 

elements were used along the crack front to approximate the r 0+1 strain singularity predicted 
from analysis. In this configuration, the nodes on the crack front are free to move independently 
while the mid-side nodes remain at the midpoints. All of the FEA were performed using 
ABAQUS. 




i 


Figure A.l. Schematic of the CC01 crack model modeled using FEA. In the present case, 
only solutions for bending (Si) were determined, the tensile load (So) was set to 
zero. 

The problem of fully plastic bending in a plate presented some challenges in the 
development of appropriate finite element models suitable for evaluating the plastic component 
of J. The first problem encountered was the formation of poorly conditioned deformed elements 
at high load levels. In fne FEA analysis performed to derive the J solutions reported in 


Appendix K, the condition for convergence of the fully plastic hj values appearing in the EPRI 
formulation for J (see Appendix K) is found by iterating the load value until the elastic J value, 
J e , is small compared to the total J value, J , along with proper (n+1) power dependence of J p on 
the load value. In the case of bending, this convergence condition proved difficult to attain. 
Near the crack front, the elements are generally small compared to the specimen dimensions, 
approximately lO^a for this study. Under the conditions of large plastic strains and high load 
values required for the condition J/J e » 1 to obtain, the deformation of the crack tip elements 
resulted in poorly conditioned elements. It was observed that the remote strain for 
configurations in which with ratio of J/Je = 100 exceeded 100%. This problem was overcome by 
performing a convergence study of the fully plastic hj value as a function of load based on J p . In 
this analysis, a finite element model was evaluated for both elastic-plastic and elastic material 
properties and J p evaluated as the J-J e . 

The second problem encountered was buckling under the applied load value. Initially, 
the bending moment was applied to the finite element model as a distributed stress along the top 
surface, using the *DLOAD user subroutine in ABAQUS. Analysis of the deformed shape of 
the finite element model showed very small crack tip opening displacement (CTOD) that varied 
little with increasing bending load. This type of deformation behavior is indicative of buckling 
in the finite element model. This problem was overcome by changing the boundary conditions 
from applied stress to applied displacement where the z-displacement was prescribed along the 
top surface of the model and the bending moment was calculated using the nodal force obtained 
from the analysis. 

The accuracy of the finite element models was examined through comparisons with 
published research. In the review of the published research, it was noted that there is little 
agreement on the fully plastic results between different authors. One of the problems noted in 
many of the publications was a lack of information regarding boundary conditions and 
convergence criteria for the fully plastic analyses. 

In the first stage of validating the FEA, the finite element models used in the current 
study were evaluated under purely linear-elastic material properties so that the results could be 
compared with the benchmark solutions of Newman & Raju. Excellent agreement was obtained 
between the current finite element model solutions and the equivalent bend solutions of Newman 
& Raju. 

In the second stage of the validation, it was hoped to compare calculated elastic-plastic J 
results with similar solutions obtained from the open literature. However, a literature review did 
not yield any published J results for the fully plastic comer crack in a plate subjected to bending. 
The closest published results found to a comer crack in bending were those for a surface crack in 
a plate under bending reported by Yagawa et al. in Three-dimensional Fully Plastic Solutions for 
Semi-elliptical Surface Cracks (Int. J of Press. Ves. and Piping, Vol. 53, pp. 457-510). As is 
evident from Figure 2, adding a symmetry boundary condition to the face with a normal in the 
positive x direction can create a surface crack model. In order to compare the present finite 
element modeling with Yagawa et al., finite element models for the surface crack in a plate 
geometry where created for two crack geometries: aJc = 0.2 and a/c = 1.0, both with aJt - 0.5. 
Significant variation was noted between the solutions of Yagawa et al. and the current finite 



element model, especially for the elongated (a/c = 0.2) crack configuration. In this case, the 
values of the crack tip parameters/; (deep crack tip, a) and/2 (surface crack tip, c) calculated by 
Yagawa et al differed by 100% and 10%, respectively, from the values generated in this study. 

Several different finite element models for the elongated crack configuration were 
created to verify the current solution. In addition to solution verification, the finite element study 
was also used to investigate the influence of mesh density and boundary conditions on the 
solution. A boundary condition of particular concern was the top surface. Yagawa et al. state 
“...axial nodal displacements along the top surface were constrained to deform linearly along the 
top surface so that it remains plane during deformation.” In the current analysis, the axial nodal 
displacements were constrained to be linear along the thickness yet the results for fully plastic ft; 
convergence yielded non-planar top surface deformations. This was caused by warping, a 
phenomenon that can accompany bending deformation and is more pronounced in those models 
with small thickness to width ratios. The results of this finite element study showed the mesh 
density in the z-direction has a small influence on /; (approximately 10%), but a much larger 
influence on fi (approximately 25%). The increased influence on / 2 is directly attributable to 
warping. 

In another set of calculations, the finite element models were constrained to reduce the 
amount of warping. The results for these cases showed an approximate 50% reduction in the 
calculated value of/;. As a result, the values of /; and /> calculated under constrained warping 
conditions now showed acceptable agreement with the results of Yagawa et al. Therefore, the 
difference between the current analysis results and those of Yagawa was demonstrated to be due 
to warping. In addition, the agreement with the results of Yagawa et al under similar boundary 
conditions validated the finite element modeling employed for the surface crack and hence, by 
implication, also the comer crack, since the surface and comer crack models only differed 
through applied boundary conditions. 

I 

The finite element models employed for CC01 geometry calculations allowed for the 
natural deformation of the specimen to occur under load controlled bending. Thus, it was not 
considered necessary to inhibit warping in these models. 

A review of FEA based J solutions in the literature revealed that there is no consistent or 
well defined method employed to define fully plastic J behavior. In the present case, the FEA 
model solutions were considered converged when successive values of ft; were within 1% for a 
constant displacement step of 0.0625 units (1.563% nominal strain), where 


ft,= 


cca 0 e o t | 


^ B+1 

kKj 


,J = J e + J f 


(Al.l) 


and M* is given by equation (7). At this point in the computations, J p was significantly larger 
than J e , and near fully plastic attained. As the loading on the finite element model produced 



nominal strains exceeding 25%, a divergence in the hj value were observed in some models. 
This divergence can be attributed to poorly conditioned deformed elements. Therefore, the finite 
element models were analyzed for convergence between 6.25% and 12.5% nominal strain. 

After the convergence load was determined, a second FEA was performed to determine 
J e at this load value. The value of hj was then calculated as 




J-J e 



(A1.2) 


The values of hj as a function of elliptical angle are shown graphically in Figures A1.2 through 
A1.7. 
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Elliptical Angle (degrees) 

Figure A1.2. Variation of hj with elliptical angle for CC01 subjected to bending, a/t=0.2, 
a/c=0.2. 




Figure A1.3. Variation of hj with elliptical angle for CC01 subjected to bending, a/t=0.2, 
a/c=0.6. 
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Elliptical Angle (degrees) 

Figure A1.6. Variation of hi with elliptical angle for CC01 subjected to bending, a/t=0.5, 
a/c=0.6. 




Elliptical Angle (degrees) 

Figure A1.7. Variation of hj with elliptical angle for CC01 subjected to bending, a/t=0.5, 
a/c=1.0. 



FEA for EC01/EC02 (Tension) Models 


In the case of the EC01/EC02 finite element model, the crack front is contained within 
the plate and there is no intersection of the crack front with a free surface. Figure A1.8 is a 
schematic of the embedded crack geometry. 



Figure A1.8. EC02 crack model. 

Symmetry allowed the geometry to be reduced to one eighth its size in the finite element 
modeling. In the FEA modeling, the ratios b/c=4 and c/h=0.25 were held constant for all the 
analyses, where b=w/2. As for the CC01 modeling, 20-noded brick elements with reduced 
integration were used with a focused ring of elements around the crack front and crack tip 

n_ 

elements that approximated the r n+1 strain singularity at the tip. 

After the convergence load was determined following similar procedures to those for the CC01 
modeling, a second FEA was performed to determine J e at this load value. The value of hj was 
then calculated as 


K =' 


J-J. 




f P y +1 


(A.3) 


The computed values of hj as a function of elliptical angle are shown graphically in Figures A1.9 
through A 1.14. 



Embedded Crack (a/c = 0.2 , a/t = 0.2) 
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Figure A1.9. Variation of hi with elliptical angle for EC01 subjected to tension, a/t=0.2, a/c=0.2. 
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Embedded Crack (a/c = 0.2, a/t = 0.5) 



Elliptical Angle (degrees) 


Figure A1.10. Variation of hj with elliptical angle for EC01 subjected to tension, a/t=0.5, 
a/c=0.2. 




Embedded Crack (a/c = 0.6, a/t = 0.2) 
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Figure Al.ll. Variation of hj with elliptical angle for EC01 subjected to tension, a/t=0.2, 
a/c=0.6. 



Embedded Crack (a/c = 0.6, a/t = 0.5) 



Figure A1.12. Variation of h] with elliptical angle for EC01 subjected to tension, a/t=0.5, 
a/c=0.6. 



Embedded Crack (a/c = 1 .0, a/t = 0.2) 



Figure A1.13. Variation of hj with elliptical angle for EC01 subjected to tension, a/t=0.2, 
a/c=1.0. 
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Embedded Crack (a/c = 1 .0, a/t = 0.5) 



Figure A1.14. Variation of h] with elliptical angle for EC01 subjected to tension, a/t=0.5, 
a/c=1.0. 




APPENDIX 2: EXAMPLE INPUT DATA 


This Appendix contains ten tables listing the data necessary to interactively input data to 
create the ten example files, Examl.inp through ExamlO.inp, contained on the distribution CD. 
The corresponding output files, Examl.out through ExamlO.out, respectively, are also contained 
on the distribution CD. Hard copies of the ten input and output files are given in Appendix 3. 


EXAMPLE 1: ELASTIC-PLASTIC J COMPUTATION 


Parameter 


Type of 
Session 


Input File 
Name 


Output File 
Name 


Parameter 


Crack Type 


Thickness 


Diameter 


Crack 

location 


Parameter 


Elastic 

Modulus 


Poisson’s 

Ratio 


Alpha 


SigmaO 


N 


Yield Stress 


Ultimate 

Stress 


Value 


5 


1 


Examplel.inp 


Example l.out 


Value 


SC 



Session Data 


Unit Description 


elastic-plastic J computation 


interactive input while creating a batch 
file 


batch file to be created 


output file for printed results 


U.S. customary units 


Crack Geomet 


Unit Description 


surface crack 


axial crack in cylinder subject to 
arbitrary hoop stressin 


inches cylinder thickness 


inches outer diameter of cylinder 


internal crack 


Material: Tensile Properties 


Value Unit Description 


30000 ksi Young’s modulus 


elastic Poisson ratio 


coefficient in Ramberg-Osgood 
equation 


yield stress in Ramberg-Osgood 
equation 


exponent in Ramberg-Osgood equation 


material 0.2% yield stress 


ultimate strength 




















































Parameter 


Loading 

Condition 


Number of 
Stresses 


Specify Unit 
Pressure? 


Pressure 



X/t 

0 

Stress 

80 

X/t 

0.25 

Stress 

60 

X/t 

0.5 

Stress 

30 

X/t 

0.75 

Stress 

0 

X/t 

1 

Stress 

-40 

X/t 

-1 


Parameter 


Print 


Master Menu? 


Quit? 


O 


Value 


1 



Applied Loadin 


Unit Description 


primary and secondary Load 


number of stress distributions to be 
specified 


Primary Load 


primary load will be specified as 
internal pressure 


ksi I internal pressure 


Secondary Load 


normalized distance 


stress value 


normalized distance 


stress value 


normalized distance 


stress value 


normalized distance 


stress value 


normalized distance 


stress value 


End Input 


Crack Sizes 


Unit Description 


manually specify crack sizes 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


crack depth 


half crack surface length 


crack depth 


half crack surface length 


crack depth 


half crack surface length 


crack depth 


half crack surface length 


crack depth 


half crack surface len 


crack depth 


half crack surface length 


End Input 


> Data 


rint results to output file 


return to master menu 


uit and save batch files 


Terminate Session 














































































EXAMPLE 2: ELASTIC-PLASTIC FAILURE ANALYSIS: CRITICAL CRACK 


Session Data 


Parameter 

Value 

Unit 

Description 

■3BSHI 

6 


elastic-plastic failure analysis 

Type of 
Session 

1 


interactive input while creating a batch 
file 

Input File 
Name 

Example2.inp 


batch file to be created 

Output File 
Name 

Example2.out 


output file for printed results 

Type of 
units 

1 


U.S. customary units 

Crack Geometry 

Parameter 

Value 

Unit 

Description 


SC 


Surface crack 

Crack Type 

2 


crack in finite width plate subject to 
arbitrary stressing 

Thickness 

1 

inches 

plate thickness 

Width 

5 

inches 

plate width 

' 

Material: Tensile Properties 

Parameter 

Value 

Unit 

Description 

Elastic 

Modulus 

30000 

ksi 

Young’s modulus 

Poisson’s 

Ratio 

0.3 


elastic Poisson ratio 

Alpha 

1 


coefficient in Ramberg-psgood 
equation 

SigmaO 

100 

ksi 

yield stress in Ramberg-Osgood 
equation 

n 

10 


exponent in Ramberg-Osgood Equation 

Yield Stress 

■BH 

ksi 

material 0.2% Yield Stress 

Ultimate 

Stress 


ksi 

ultimate strength 

Crack Shape 

Parameter 

Value 

Unit 

Description 

Constant 

Aspect 

Ratio? 

1 


crack has constant aspect ratio as 
opposed to constant surface length 


0.35 


value of a/c 




























































Analysis Type 


Parameter 

Value 

Unit 

Description 

Critical 

Crack 

1 


perform critical crack calculation as 
opposed to critical load 

Ductile 

2 


Perform ductile analysis as opposed to 
brittle analysis 


Parameter 

Value 

Unit 

Description 

Toughness, 

Jmat 

0.25 

ksi-inch 

Toughness expressed in terms of the J- 
integral 

Resistance 

Curve 

1 


The J-R curve is expressed as a 
quadratic form in the tear length 

DjO 

0.245 

ksi -inch 

First coefficient of quadratic equation 
for J R 

m 

30 

ksi 

Coefficient of linear J-R term 

Dj2 

-50 


Coefficient of quadratic J-R term 

dflniax 

0.3 

inches 

Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 


Parameter 

Value 

Unit 

Description 

Interactively 
input load? 

1 


The loads will be specified manually 

Loading 

Condition 

2 


Primary and Secondary Load 

Number of 
Stresses 

1 


Number of Stress Distributions to be 
Specified 


Primary Load 


x/t 

0 


Normalized Distance 

Stress 

120 

ksi 

Stress Value 

x/t 

1 


Normalized Distance 

Stress 

120 

ksi 

Stress Value 

x/t 

-1 


End Input 


Secondary Load 


x/t 


0 


Normalized Distance 


Stress 


100 


ksi 


Stress Value 


x/t 


Normalized Distance 


Stress 


50 


-1 


ksi 


Stress Value 


End Input 


x/t 



















































































EXAMPLE 3: ELASTIC-PLASTIC FAILURE ANALYSIS: CRITICAL LOAD 


Parameter 

Value 

HBH! 

6 

Type of 
Session 

1 

Input File 
Name 

Example3.inp 

Output File 
Name 

Example3.out 

Type of 
units 

2 


Parameter 

Value 


EC 

Crack Type 

2 

Thickness 

0.05 

Width 

0.15 


Parameter 

Value 

Elastic 

Modulus 

210000 

Poisson’s 

Ratio 

0.3 

Alpha 

1.5 

SigmaO 

400 

n 

15 

Yield Stress 

400 

Ultimate 

Stress 

600 


Parameter Value 

Constant 1 

Aspect 

Ratio? 




Session Data 



Description 

elastic-plastic failure analysis 

interactive input while creating a batch 

file 

batch file to be created 


output file for printed results 


Metric units 


Crack Geomet 




Description 

embedded crack 


crack in finite width plate subject to 
arbitrary stressing 



meters 

plate thickness 


meters 

plate width 


Material: Tensile Properties 


VnU 

Description 


MPa 

Young’s modulus 



elastic Poisson ratio 


coefficient in Ramberg-Osgood 

equation 

MPa yield stress in Ramberg-Osgood 




equation 



exponent in Ramberg-Osgood Equation 


IHESflH 

material 0.2% Yield Stress 


MPa 

ultimate strength 


Crack Shape 



Unit 

Description 



crack has constant aspect ratio as 
opposed to constant surface length 


0.25 


value of a/c 





























































Analysis Type 


Parameter 

Value 

Unit 

Description 

Critical 

Loads 

2 


perform critical load calculation as 
opposed to critical crack 

Ductile 



2 


Perform ductile analysis as opposed to 
brittle analysis 

Material Toughness Properties 

Parameter 

Value 

Unit 

Description 

Toughness, 

Jmat 

0.05 

MPa-meter 

Toughness expressed in terms of the J- 
integral 

Resistance 

Curve 

2 


The J-R curve is expressed as the tear 
length raised to a power 

Dji 

5.0 

MPa-meter 

coefficient of power law equation for Jr 

Dj2 

0.5 

MPa 

exponent of power law for J-R term 

d^max 

0.004 

meters 

Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 

Loading Condition 

Parameter 

Value 

Unit 

Description 

Loading 

Condition 

0 


Primary load only 

Number of 
Stresses 

1 


Number of Stress Distributions to be 
Specified 


x/t 

0 

Stress 

400 

x/t 

1 

Stress 

400 

x/t 

-1 



Normalized Distance 

Stress Value 

Normalized Distance 
Stress Value 



End Input 



Crack Sizes 


Value 

Unit 

Description 

1 


manually specify crack sizes 

0.003 

meters 

crack depth 

0.004 

meters 

half crack surface length 

0.005 

meters 

crack depth 

-1 


End input 


Session Data 


sto 



End Input Data 












































































EXAMPLE 4: ELASTIC-PLASTIC FATIGUE LIFE ANALYSIS 


Session Data 


Parameter 

Value 

Unit 

Description 


1 


elastic-plastic failure analysis 

Type of 
Session 

1 


interactive input while creating a batch 
file 

Input File 
Name 

Example4.inp 


batch file to be created 

Output File 
Name 

Example4.out 


output file for printed results 

Type of 
units 

1 


U.S. customary units 

Crack Geometry 

Parameter 

Value 

Unit 

Description 


CC 


Comer crack 


1 


crack in finite width plate 

Tension 

1 


Tension as opposed to bending 

Thickness 

2. 

inches 

plate thickness 

Width 

2 

inches 

plate width 


Material: Tensile Properties 

Parameter 

Value 

Unit 

Description 

Elastic 

Modulus 

30000 

ksi 

Young’s modulus 

Poisson’s 

Ratio 

0.3 


elastic Poisson ratio 

Alpha 

1 


coefficient in Ramberg-Osgood 
equation 

SigmaO 

100 

ksi 

yield stress in Ramberg-Osgood 
equation 

n 

10 


exponent in Ramberg-Osgood Equation 

Yield Stress 

100 

ksi 

material 0.2% Yield Stress 

Ultimate 

300 

ksi 

ultimate strength 


Stress 


Service Load Spectrum (Schedule) 


Block Case Definition: Case # 1 : Maximum Load Data 


Maximum 
Load Type 


Primary and Secondary load 


Primary Maximum Load 


ksi 


Primary tensile stress at maximum load 


Tensile 

stress 


80 




























































Secondary Maximum Load 


x/t 

0 


Normalized Distance 

Stress 

50 

ksi 

Stress Value 

x/t 

1 


Normalized Distance 

Stress 

50 

ksi 

Stress Value 

x/t 

-1 


End Input 


Block Case Definition: Case #h Minimum Load Data 


Minimum 
Load Type 


Primary and Secondary load 


Primary Minimum Load 


Tensile 

stress 


0 


ksi 


Primary tensile stress at maximum load 


Secondary Minimum Load 


x/t 

0 


Normalized Distance 

Stress 

50 

ksi 

Stress Value 

x/t 

1 


Normalized Distance 

Stress 

50 

ksi 

Stress Value 

x/t 

-1 


End Input 


Cycles 


Number of 
cycles 


Number of times maximum and 
minimum loads for Case # 1 are applied 


Block Case Definition: Case #2: Maximum Load 


Maximum 
Load Type 


Primary and Secondary load 


Primary Maximum Load 


Tensile 


50 


ksi 


Primary tensile stress at maximum load 


stress | 


Secondary Maximum Load 

x/t 

0 


Normalized Distance 

Stress 

50 

ksi 

Stress Value 

x/t 

1 


Normalized Distance 

Stress 

50 

ksi 

Stress Value 

x/t 

-1 


End Input 

Block Case Definition: Case #2: Maximum Load 

Minimum 
Load Type 

2 


Primary and Secondary load 

Primary Minimum Load 

Tensile 0 

ksi 

Primary tensile stress at maximum load 


Tensile 

stress 


0 































































x/t 

0 

Stress 

50 

x/t 

1 

Stress 

50 

x/t 

-1 




Block Case 

# 


Block Case 
# 


Parameter 


Toughness, 

Jmat 


Fatigue 

coefficient 



Uo 

1 

Alp (a-tip) 

1 

Alp (c-tip) 

1 


Secondary Minimum Load 


Normalized Distance 


Stress Value 


Normalized Distance 


Stress Value 


End Input 


Cycles 


Number of times maximum and 
minimum loads for Case # 2 are applied 


Block Case Definition: Case #3: Maximum Load 


Terminate Block Case Input 


Load Spectrum (Schedule . 


1 Block Case number or ID 


Number of times Block Case 1 is 
applied 


Block Case number or ID 


Number of times Block Case 2 is 
applied 


End load spectrum input 


Material Properties 


Unit Description 


ksi-inch Toughness expressed in ( terms of the J- 

integral ; 


Chosen so crack Coefficient in Paris equation 
growth rate is in 
inches/cycle 


Exponent in Paris equation 


Crack closure term appropriate to Paris 
equation test data 


Constraint factor for a-tip (used in crack 
closure evaluation) 


Constraint factor for c-tip (used in crack 
closure evaluation) 


Initial Crack Size 


Description 



Value 


0.5 



Parameter 

Value 

Unit 

a 

0.05 

inches 

c 

0.1 

inches 



























































Maximum 
number of 
schedules 

100 

Print 

interval 

10 


Print 

p 

Master 

Menu? 

0 

■BsSHiSi 

i 

■EKBI 

0 


Schedule and Print Data 


The fatigue calculations will terminate 
after the load schedule has been applied 
this number of times if failure has not 
occurred first 


Results will be printed to the output file 
after the schedule has been applied this 
number of times 


Post Analysis Data 


Print Results to Output File 


Return to Master Menu 


Quit and Save Batch Files 


Terminate Session 
















EXAMPLE 5: PROOF TEST PROCEDURE: PRE-PROOF TEST ANALYSIS: 

CRITICAL CRACK 


Session Data 


Parameter 

Value 

Unit 

Description 

■333391 

8 


elastic-plastic failure analysis 

Type of 
Session 

1 


interactive input while creating a batch 
file 

Input File 
Name 

Example5.inp 


batch file to be created 

Output File 
Name 

Example5.out 


output file for printed results 

Type of 
units 

1 


U.S. customary units 

Proof Test 
Analysis 

1X2? 

1 


Perform Pre-Proof Test Safe Life 
Analysis as opposed to Proof Test 
Analysis 

Pre-Proof 

Analysis 

Tm 

1 


Perform Critical Flaw Size calculation 
as opposed to Fatigue Life calculation. 


Crack Geometry 


Parameter 

Value 

Unit 

Description 


SC 


Surface crack 

Crack Type 

2 


crack in finite width plate subject to 
arbitrary stressing 

Thickness 

1 

inches 

plate thickness 

Width 

5 

inches 

plate width 


Material: Tensile Properties 


Parameter 

Value 

Unit 

Description 

Elastic 

Modulus 

30000 

ksi 

Young’s modulus 

Poisson’s 

Ratio 

0.3 


elastic Poisson ratio 

Alpha 

1 


coefficient in Ramberg-Osgood 
equation 

SigmaO 

100 

ksi 

yield stress in Ramberg-Osgood 
equation 

n 

10 


exponent in Ramberg-Osgood Equation 

Yield Stress 

100 

ksi 

material 0.2% Yield Stress 

Ultimate 

Stress 

300 

ksi 

ultimate strength 


























































Crack Shape 


Parameter 

Value 

Unit 

Description 

Constant 

Aspect 

Ratio? 

1 


crack has constant aspect ratio as 
opposed to constant surface length 


0.6 


value of a/c 

Analysis Type 

Parameter 

Value 

Unit 

Description 

Ductile 

2 


Perform ductile analysis as opposed to 
brittle analysis 

N 

laterial Toughness Properties 

Parameter 

Value 

Unit 

Description 

Toughness, 

Jmat 

0.25 

ksi-inch 

Toughness expressed in terms of the J- 
integral 

Resistance 

Curve 

1 


The J-R curve is expressed as a 
quadratic form in the tear length 

DjO 

0.245 

ksi-inch 

First coefficient of quadratic equation 
for J R 

Djl 

30 

ksi 

Coefficient of linear J-R term 

Dj2 

-50 

Ktxsaai 

Coefficient of quadratic J-R term 

d&max 

0.3 

inches 

Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 

Applied Service Loading 

Parameter 

Value 

Unit 

Description 

Manually 
input load? 

1 


The loads will be specified manually 

Loading 

Condition 

2 


Primary and Secondary Load 

Primary Service Load 

x/t 

0 


Normalized Distance 

Stress 

no 

ksi 

Stress Value 

x/t 

l 


Normalized Distance 

Stress 

120 

ksi 

Stress Value 

x/t 

-1 


End Input 


Secondary Service Load 


x/t 

0 

Stress 

100 

x/t 

1 

Stress 

50 

x/t 

-1 



Normalized Distance 

Stress Value 

Normalized Distance 

Stress Value 

End Input 





















































































EXAMPLE 6: PROOF TEST PROCEDURE: FLAW SCREENING 


Parameter 


Type of 
Session 


Input File 
Name 


Output File 
Name 


Type of 
units 



Example6.inp 


Example6.out 



Flaw 

Screening 


Parameter 



Crack Type 


Thickness 


Diameter 


External 


Parameter 


Elastic 

Modulus 


Poisson’s 

Ratio 


Alpha 


SigmaO 


Yield Stress 


Ultimate 

Stress 


Value 


SC 



0.084 


1.183 


Value 


210000 


Session Data 


Unit Description 


elastic-plastic failure analysis 


interactive input while creating a batch 
file 


batch file to be created 


output file for printed results 


Metric units 


Perform Proof Test Analysis as opposed 
to Pre-Proof Test Safe Life Analysis 


Perform Flaw Screening Analysis as 
opposed to Proof Load Analysis or Final 
Crack Size Analysis 


Crack Geometr 


Description 


Surface crack 


crack in cylinder subject to arbitrary 
hoop stressin 


cylinder thickness 


Outer diamater of cylinder 


Crack is on the outside of the cylinder 


Material: Tensile Properties 


Unit Description 


MPa Young’s modulus 


elastic Poisson ratio 


coefficient in Ramberg-Osgood 
equation 


yield stress in Ramberg-Osgood 
equation 


exponent in Ramberg-Osgood Equation 


material 0.2% Yield Stress 


ultimate strength 


meters 


meters 



























































Crack Shape ! 

Parameter 

Value 

Unit 

Description 


0.61 


value of a/c 

Analysis Type 

Parameter 

Value 

Unit 

Description 

Ductile 

2 


Perform ductile analysis as opposed to 
brittle analysis 

Material Toughness Properties | 

Parameter 

Value 

Unit 

Description 

Toughness, 

Jmat 

0.433 

MPa-meter 

Toughness expressed in terms of the J- 
integral 

Resistance 

Curve 

1 

' 

The J-R curve is expressed as a 
quadratic form in the tear length 

DjO 

0.361 

MPa-meter 

First coefficient of quadratic equation 
for J R 

5iL 

503.3 

MPa 

Coefficient of linear J-R term 

Dj2 

-2325 


Coefficient of quadratic J-R term 

dSrnax 

0.005 

meter 

Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 

Applied Proof Loading 

Parameter 

Value 

Unit 

Description 

Manually 
input load? 

1 


The loads will be specified manually 

Loading 

Condition 

1 


Primary Load ^ 

| Primary Load 

Specify Unit 
Pressure? 

Y 


primary load will be specified as 
internal pressure 

Pressure 

80.45 

MPa 

internal pressure 

Proof Load 
Factor 

1 


The proof load is scaled by this factor 

! Post Analysis Data 

Print 

P 


Print Results to Output File 


0 


Return to Master Menu 

| Quit? 

1 


Quit and Save Batch Files 


Terminate Session 























































EXAMPLE 7: PROOF TEST PROCEDURE: PROOF LOAD 


Parameter 


Type of 
Session 


Input File 
Name 


Output File 
Name 


Type of 
units 


Value 


8 


1 


Example7.inp 


Example7.out 



Parameter 


wassssm 


Crack Type 
Thickness 


Width 


Parameter 


Elastic 

Modulus 


Poisson’s 

Ratio 


Alpha 


SigmaO 


Yield Stress 


Ultimate 

Stress 


Parameter 


Value 


EC 


Session Data 


Unit Description 


elastic-plastic failure analysis 


interactive input while creating a batch 
file 


batch file to be created 


output file for printed results 


U.S. Customary units 


Perform Proof Test Analysis as opposed 
to Pre-Proof Test Safe Life Analysis 


Perform Proof Load Analysis as 
opposed to Flaw Screening Analysis or 
Final Crack Size Analysis 


Crack Geomet 


Unit Description 


Embedded crack 



inches 


meters 


Value 


30000 


late thickness 


Plate width 


Material: Tensile Properties 


Unit Description 


ksi Young’s modulus 


elastic Poisson ratio 



Value 


0.25 


coefficient in Ramberg-Osgood 
equation 


yield stress in Ramberg-Osgood 
equation 


exponent in Ramberg-Osgood 


material 0.2% Yield Stress 


ultimate strength 


Crack Shape 


Unit Description 


value of a/c 


uation 































































Parameter 

Value 

Ductile 

2 


Analysis Type 


Unit 


Parameter 


Toughness, 

Jmat 


Resistance 

Curve 


Value 


0.2 


Description 



Parameter 

Value 

Loading 

Condition 

1 


X/t 

0 

stress 

60 

X/t 

0.2 

stress 

55 

X/t 

0.4 

stress 

52 

X/t 

0.6 

stress 

50 

X/t 

0.8 

stress 

49 

X/t 

1 

stress 

47 

X/t 

-1 


X/t 

0 

stress 

90 

X/t 

1 

stress 

0 

X/t 

-1 


Perform ductile analysis as opposed to 
brittle analysis 


Material Toughness Properties 


Unit Description 


ksi-inch Toughness expressed in terms of the J- 
integral 


The J-R curve is expressed as the tear 
length raised to a power 


Coefficient of linear J-R term 


Coefficient of quadratic J-R term 


inches Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 


Applied Proof Loadin 


Description 


Primary and Secondary Loads present 
during proof testing 


Primary Proof Load 


Normalized distance 


Stress value 


Normalized distance 


Stress value 


Normalized distance 


Stress value 


Normalized distance 


Stress value 


Normalized distance 


Stress value 


Normalized distance 


Stress value 


End input 


Secondary Proof Load 


Normalized distance 


ksi Stress value 


Normalized distance 


ksi Stress value 


End input 
























































































































EXAMPLE 8: PROOF TEST PROCEDURE: FINAL CRACK SIZE 


Session Data 


Parameter 

Value 

Unit 

Description 


8 


elastic-plastic failure analysis 

Type of 
Session 

1 


interactive input while creating a batch 
file 

Input File 
Name 

Example8.inp 


batch file to be created 

Output File 
Name 

Example8.out 


output file for printed results 

Type of 
units 

1 


U.S. Customaryc units 

Proof Test 
Analysis 
Type . 

3 


Perform Proof Test Analysis as opposed 
to Pre-Proof Test Safe Life Analysis 

Final Crack 
Size 

3" 


Perform Final Crack Size Analysis as 
opposed to Proof Load Analysis or Flaw 
Screening Analysis 

Crack Geometry 

Parameter 

Value 

Unit 

Description 

Emssm 

CC 


Surface crack 


1 


Comer crack in plate 

Bending 

2 


Bending in the thickness direction as 
opposed to tension 

Thickness 

2 

inches 

plate thickness 

Width 

2 

inches 

Plate width 


Material: Tensile Properties | 

Parameter 

Value 

Unit 

Description 

Elastic 

30000 

ksi 

Young’s modulus 


Modulus 

Poisson’s 

Ratio 


0.3 


elastic Poisson ratio 


Alpha 


1.2 


coefficient in Ramberg-Osgood 
equation 


SigmaO 


150 


ksi 


yield stress in Ramberg-Osgood 
equation 


n 


20 


exponent in Ramberg-Osgood Equation 


Yield Stress 


150 


MPa 


material 0.2% Yield Stress 


Ultimate 

Stress 


200 


MPa 


ultimate strength 






















































Crack Shape 

Parameter 

Value 

Unit 

Description 


0.5 

inches 

Initial crack depth at the a-tip 


1 


Initial aspect ratio, a/c 

Proof Load 

iaia 

0 


No secondary stress is present during 
the proof test 


220 

ksi 

Primary proof load 

Material Toughness Properties 

Parameter 

Value 

Unit 

Description 

Resistance 

Curve 

1 


The J-R curve is expressed as a 
quadratic form in the tear length 

DjO 

0.145 

Ksi-inch 

First coefficient of quadratic equation 
for J R 


Djl 

30 

ksi 

Coefficient of linear J-R term 

Dj2 

-50 

besesssbi 

Coefficient of quadratic J-R term 

d^max 

0.3 

inch 

Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 

Toughness, 

Jmat 

0.15 

Ksi-inch 

Toughness expressed in terms of the J- 
integral 


Post Analysis Data 


Print 


Print Results to Output File 


Master 

Menu? 


0 


Return to Master Menu 


Quit? 


Quit and Save Batch Files 


Option 


0 


Terminate Session 










































EXAMPLE 9: ELASTIC-PLASTIC TEAR-FATIGUE LIFE ANALYSIS 


Parameter 



Type of 
Session 

1 

Input File 
Name 

Example9.inp 

Output File 
Name 

Example9.out 

Type of 
units 

1 


Parameter 


Value 


SC 


Session Data 


Unit Description 


elastic-plastic failure analysis 


interactive input while creating a batch 
file 


batch file to be created 


output file for printed results 


U.S. customary units 


Crack Geometr 


Unit 


Surface crack 




4 


Axial crack in cylinder 

Thickness 

0.5 

inches 

Thickness of cylinder 

Diameter 

40 

inches 

Outer diameter of cyklinder 

Crack 

location 

e 


External crack 

1 

Material: Tensile Properties 

Parameter 

Value 

Unit 

Description 

Elastic 

Modulus 

30000 

ksi 

Young’s modulus 

Poisson’s 

Ratio 

0.3 


elastic Poisson ratio 

Alpha 

1 


coefficient in Ramberg-Osgood 
equation 

SigmaO 

100 

ksi 

yield stress in Ramberg-Osgood 
equation 

N 

10 


exponent in Ramberg-Osgood Equation 

Yield Stress 

100 

ksi 

material 0.2% Yield Stress 

Ultimate 

Stress 

150 

ksi 

ultimate strength 

Crack Size 

Crack depth 

0.25 


Depth of surface crack 


0.5 


a/c value 































































Service Load Histo 


Number of service cycles 


Primary and Secondary load 


Primary Maximum Load 

(Note : Stress distributions are specified with respect to the inner surface of the cylinder 

for external cracks.) 


Primary stress distribution is not 



identical to that due to internal pressure 


Normalized Distance 


Stress Value 


Normalized Distance 


Stress Value 


End Input 


Secondary maximum load 

(Note: Stress distributions are specified with respect to the inner surface of the cylinder 

'or external cracks.) 


Normalized Distance 


ksi Stress Value 


Normalized Distance 


ksi Stress Value 


End Input 


Primary Minimum Load 

(Note: Stress distributions are specified with respect to the irmer surface of the cylinder 

for external cracks.) 


Primary stress distribution is not 


x/t 


Stress 


x/t 


Stress 


x/t 


, x/t 


Stress 


x/t 


Stress 



identical to that due to internal pressure 


x/t 0 Normalized distance 


Stress 0 ksi Stress Value 


x/t 1 Normalized distance 


Stress 0 ksi Stress Value 


x/t -1 End Input 


Secondary minimum load 

(Note: Stress distributions are specified with respect to the inner surface of the cylinder 

for external cracks.) 


x/t 0 Normalized Distance 


Stress -100 ksi Stress Value 


x/t 1 Normalized distance 


stress 100 ksi Stress value 


x/t -1 End Input 



































































Material Properties | 

Parameter 

Value 

Unit 

Description 

Resistance 

Curve 

1 


The J-R curve is expressed as a 
quadratic form in the tear length 

DjO 

0.2 

Ksi-inch 

First coefficient of quadratic equation 
for J R 

Djl 

30 

ksi 

Coefficient of linear J-R term 

m 

0 

WMMEEMM 

Coefficient of quadratic J-R term 

d^max 

0.2 

inch 

Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 

Toughness, 

Jmat 

0.21 

Ksi-inch 

Toughness expressed in terms of the J- 
integral 

U 0 

1 


Crack closure term appropriate to Paris 
equation test data 

Alp (a-tip) 

1 


Constraint factor for a-tip (used in crack 
closure evaluation) 

Alp (c-tip) 

1 


Constraint factor for c-tip (used in crack 
closure evaluation) 

Fatigue 

coefficient 

W™ 

Chosen so crack 
growth rate is in 
inches/cycle 

Coefficient in Paris equation 

Fatigue 

exponent 

4 


Exponent in Paris equation 

I Post Analysis Data 

Print 

P 


Print Results to Output File 

Master 

Menu? 

0 


Return to Master Menu* 

i 

Quit? 

1 


Quit and Save Batch Files 


Op 



Terminate Session 












































EXAMPLE 10: MULTI CYCLE PROOF TEST ANALYSIS 


Session Data 


Parameter 

Value 

Unit 

Description 

wsmm 

10 


elastic-plastic failure analysis 

Type of 
Session 

1 


interactive input while creating a batch 
file 

Input File 
Name 

ExamplelO.inp 


batch file to be created 

Output File 
Name 

ExamplelO.out 


output file for printed results 

Type of 
units 

1 


U.S. customary units 

Crack Geometry 

Parameter 

Value 

Unit 



SC 


Surface crack 


2 


Crack in plate 

Thickness 

1 

inches 

plate thickness 

Width 

10 

inches 

plate width 


Material: Tensile Properties 

Parameter 

Value 

Unit 

Description 

Elastic 

Modulus 

30000 

ksi 

Young’s modulus 

Poisson’s 

Ratio 

0.3 


elastic Poisson ratio 

Alpha 

1 


coefficient in Ramberg-Osgood 
equation 

SigmaO 

100 

ksi 

yield stress in Ramberg-Osgood 
equation 

n 

5 


exponent in Ramberg-Osgood Equation 

Yield Stress 

100 

ksi 

material 0.2% Yield Stress 

Ultimate 

Stress 

200 

ksi 

ultimate strength 

Crack Size 

Distribution 

constant 

0.15 


Constant in the crack size exponential 
distribution function 


0.5 


a/c value 


Proof and Service Load Histories 


Number of times proof load is applied 
Number of service cycles 


Proof cycles 


Cycles 


50 



Primary load only 






























































Primary Maximum Proof Load 

x/t 

0 


Normalized Distance 

Stress 

120 

ksi 

Stress Value 

x/t 

1 


Normalized Distance 

Stress 

120 

ksi 

Stress Value 

x/t 

-1 


End Input 


Primary Minimum Proof Load 

x/t 

0 


Normalized Distance 

Stress 

0 

ksi 

Stress Value 

x/t 

1 


Normalized Distance 

Stress 

0 

ksi 

Stress Value 

x/t 

-1 


End Input 

Service Load History (Schedule) 

Number of 
load blocks 

1 


Number of load blocks that constitute 
the service load history 

cycles 

2 


Number of cycles load block 1 is 
applied 




































































Parameter 


Resistance 

Curve 


DjO 


Toughness, 

Jmat 


U 0 


Alp (a-tip) 


Alp (c-tip) 


Proof 

Fatigue 

coefficient 


Service 

Fatigue 

coefficient 


Print 


Master 

Menu? 


Quit? 


Value 


1 


Material Properties 


Unit 


Ksi-inch 


Ksi-inch 


Chosen so crack 
growth rate is in 
inches/cycle 


Description 


The J-R curve is expressed as a 
quadratic form in the tear length 


First coefficient of quadratic equation 
for J R 


Coefficient of linear J-R term 


Coefficient of quadratic J-R term 


Saturation tear length, the value of Jr is 
constant for tear lengths that exceed this 
value 


Toughness expressed in terms of the J- 
integral 


Crack closure term appropriate to Paris 
equation test data 


Constraint factor for a-tip (used in crack 
closure evaluation) 


Constraint factor for c-tip (used in crack 
closure evaluation) 


Coefficient in Paris equation applicable 
to proof test 


Exponent in Paris equation applicable to 
proof test 


Service fatigue crack growth properties 
same as proof test properties 


Chosen so crack Coefficient in Paris equation applicable 
growth rate is in to proof test 
inches/cycle 


Exponent in Paris equation applicable to 
proof test 


Post Analysis Data 


Print Results to Output File 


Return to Master Menu 


Quit and Save Batch Files 


Terminate Session 












































APPENDIX 3: INPUT AND OUTPUT FILES FOR EXAMPLES 1 THROUGH 12 


EXAMPLE 1 


INPUT FILE: Examl.inp 

Examl.out Output file name*12 

1 1=US units; 2=SI units 
sc Crack Model Type 

4 Crack Model Number 

0 . 250000E+00 T 

0.600000E+01 Outer Diameter 

i 

0.3000E+05 Elastic Young's modulus 

0.300 Poisson's ratio 

1.000 Alpha 

0 . 1000E+03 SigmaO 

10.000 n 

0 . 1000E+03 material yield stress 
0.1500E+03 material ultimate stress 
2 1: Primary, 2: Primary+Secondary 


1 # of Stress Dist 
y Interal Pressure (Y/N) 


10.00 

Internal Pressure 

0 . 0000E+00 

Non-Dimensional 

position 

0 . 8000E+02 

Stress value 


0 . 2500E+00 

Non-Dimensional 

position 

0 . 6000E+02 

Stress value 


0 . 5000E+00 

Non-Dimensional 

position 

0 . 3000E+02 

Stress value 


0 . 7500E+00 

Non-Dimensional 

position 

0 . 0000E+00 

Stress value 


0 . 1000E+01 

Non-Dimensional 

position 

- . 4000E+02 

Stress value 


- . 1000E+01 

Non-Dimensional 

position 

1 l=interactively input, : 

2=tabulate the data 

0 . 0500E+00 

a ( 1) 


0 . 0500E+00 

c ( 1) 


0 . 0500E+00 

a ( 2) 


0 . 1000E+00 

c ( 2) 


0 . 0500E+00 

a ( 3) 


0 . 1500E+00 

c ( 3) 


0 . 1000E+00 

a ( 4) 


0 . 1000E+00 

c( 4) 


0 . 1000E+00 

a ( 5) 


0 . 1500E+00 

c ( 5) 


0 . 1500E+00 

a ( 6) 


0 . 1500E+00 

c ( 6) 


- . 1000E+01 

a ( 7) 


P P(lst 

col . ) : to print 



0 1 : to resume, 0: stop 



OUTPUT FILE: Examl.out 


ELASTIC-PLASTIC J CALCULATION FOR SC04 


DATE: 18-SEP-03 TIME: 15:17:22 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = examl.inp 
Output Filename = Examl.out 

Cylinder Thickness, t = 0.2500 

Outer Diameter, D = 6.0000 

Crack Type = INTERNAL 

Material Yield Stress = 100.00 

Material Ultimate Stress = 150.00 

Data for the Nonlinear Material Behavior: 

SigO = 0 . 1000E+03 

E = 0.3000E+05 

nu = 0 . 3000E+00 

alpha = 0 . 1000E+01 

n = 0 . 1000E+02 

Internal Pressure 10.000 
PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

Stresses at 10 points 

Distance x/t Stress due to Total Stress 

(from inner wall) Int. Pressure 


0.0000 

125.2174 

125.2174 

0.1111 

123.9715 

123.9715 

0.2222 

122.7624 

122.7624 

0.3333 

121.5887 

121.5887 

0.4444 

120.4490 

120.4490 

0.5556 

119.3420 

119.3420 

0.6667 

118.2665 

118.2665 

0.7778 

117.2214 

117.2214 

0.8889 

116.2053 

116.2053 

1.0000 

115.2174 

115.2174 

SECONDARY LOAD 

DISTRIBUTION: 


Norm, x 

Stress 


0.00 

0 . 8000E+02 


0.25 

0 . 6000E+02 


0.50 

0 . 3000E+02 


0.75 

0 . 0000E+00 


1.00 

- . 4000E+02 




Model Code= SC04 


a 

J(c) 


c 


Je(a) Je(c) Jp(a) Jp(c) J(a) 


0 . 500E-01 
0 . 227E+00 
0 . 500E-01 
0 . 221E+00 
0 . 500E-01 
0 . 204E+00 
0 . 100E+00 
0 . 500E+00 
0 . 100E+00 
0 . 540E+00 
0 . 150E+00 
0 . 895E+00 


0 . 500E-01 
0.100E+00 
0 . 150E+00 
0 . 100E+00 
0.150E+00 
0 . 150E+00 


0 . 104E+00 
0 . 164E+00 
0.203E+00 
0 . 189E+00 
0 . 265E+00 
0 . 254E+00 


0 . 125E+00 
0 . 127E+00 
0.121E+00 
0.277E+00 
0.298E+00 
0.476E+00 


0 . 591E-01 
0 . 110E+00 
0 . 136E+00 
0 . 123E+00 
0 . 193E+00 
0 . 211E+00 


0 . 102E+00 
0 . 943E-01 
0 . 825E-01 
0 . 223E+00 
0 . 242E+00 
0 . 418E+00 


0 . 164E+00 
0 . 274E+00 
0 . 339E+00 
0 . 312E+00 
0 . 458E+00 
0 . 466E+00 


EXAMPLE 2 


INPUT FILE: Exam2.inp 


Exam2.out Output file name*12 

1 1=US units; 2=SI units 
sc Crack Model Type 

2 Crack Model Number 
0.1000E+01 Thickness 
0.5000E+01 Width 

0.3000E+05 Elastic Young's modulus 

0.300 Poisson's ratio 

1.000 Alpha 

0.1000E+03 SigmaO 

10.000 n 


0.1000E+03 material yield stress 
0.2000E+03 material ultimate stress 

1 1: const, asp. 2: const, length 

0.3500E+00 constant aspect ratio 

1 1: crit. crack 2: crit. load 

2 1: brittle, 2: ductile 

0.2500E+00 matl toughness 

1 1 : quad . 2 : power 

0.2450E+00 djO — quadratic 

0.3000E+02 djl — quadratic 

-.5000E+02 dj2 — quadratic 

0.3000E+00 da (max) — quadratic 

L 1= interactively input, 2= create a table 

2 1= pri., 2=pri. & sec. 

1 # of Stress Dist 


O.OOOOOOOOOOOOOOOE+OOO 

120.000000000000 

1.00000000000000 

120.000000000000 

- 1.00000000000000 


Nondim position 
Stress value 
Nondim position 
Stress value 
Nondim position 



0 . OOOOE+OO Non-Dimensional position 
0.1000E+03 Stress value 
0.1000E+01 Non-Dimensional position 
0.5000E+02 Stress value 
- . 1000E+01 Non-Dimensional position 
1.000 Load Factor # 1 

0 1 = input, 0 = stop 

P P(lst col.): to print 

0 1 : to resume, 0: stop 


OUTPUT FILE: Examlout 

ELASTIC- PLASTIC ANALYSIS FOR CRITICAL CRACK/LOAD FOR SC02 


DATE: 17-SEP-03 TIME: 16:04:43 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = Exam2 . inp 

Output Filename = Exam2.out 

Plate Thickness, t = 1.0000 

’ Width, W = 5.0000 

Material Yield Stress = 100.00 

Material Ultimate Stress = 200.00 

Data for the Nonlinear Material Behavior: 

SigO = 0 . 1000E+03 

E = 0 . 3000E+05 

nu = 0 . 3000E+00 

alpha = 0 . 1000E+01 

n = 0 . 1000E+02 

Data for the Elastic Plastic Failure Analysis 
♦DUCTILE ANALYSIS* is performed 
Ultimate Tensile Stress (Su) = 0.2000E+03 
Jmat = 0 . 2500E+00 

Kmat(c) = 0 . 8660E+02 , Kmat(a) = 0.9078E+02 

Search for * CRITICAL CRACK LENGTH* 

Constant aspect ratio - 0.3500E+00 

Fracture resistance curve (quadratic form): 

Jr = (0.2450E+00)+(0.3000E+02) *x+ (-.5000E+02) *x~2 

da(max) = 0.3000E+00 

Model Code= SC02 
PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

0.00 0 . 1200E+03 



1.00 0 . 1200E+03 

SECONDARY LOAD DISTRIBUTION: 

Norm, x Stress 

0.00 0 . 1000E+03 

1.00 0 . 5000E+02 

TABLE OUTPUT (Fixed a/c=0 . 350E+00) : 

Pri. Load a_iiiit a_crit a_inst da (tear) P/P0_init P/P0_crit 

P/P0_inst 

c_init c_crit c_inst dc(tear) 

0 . 100E+01 0 . 163E-01 0.111E+00 0.257E+00 0.147E+00 0.120E+01 0.121E+01 
0.124E+01 

0 . 467E-01 0 . 316E+00 0.416E+00 0.101E+00 


EXAMPLE 3 


INPUT FILE : Exam3,inp 

Exam3 . out Output file name*12 

2 1=US units; 2=SI units 
ec Crack Model Type 

2 Crack Model Number 

0 . 500E-01 Thickness 
0 . 150E+00 Width 

0 . 2100E+06 Elastic Young" s modulus 

0.300 Poisson"s ratio 
1.500 Alpha 
0.4000E+03 SigmaO 

15.000 n 

0.4000E+03 material yield stress 
0.6000E+03 material ultimate stress 

1 1: const, asp. 2: const, length 

0.2500E+00 constant aspect ratio 

2 1: crit. crack 2: crit. load 

2 1: brittle, 2: ductile 

0.5000E-01 matl toughness 

2 1 : quad . 2 ; power 

0.5000E+01 djl ~ power law 

0.5000E+00 dj2 -- power law 

0.4000E-02 da (max) -- power law 

0 Is with 2nd load, 0: w/o 2nd load 

1 # of Stress Dist 

0 . OOOOOOOOOOOOOOOE+OOO Nondim position 

400.000000000000 Stress value 

1.00000000000000 Nondim position 

400.000000000000 Stress value 

-1.00000000000000 Nondim position 

1 1= interactively input, 2=create a table 

0 . 3000E-02 a ( 1) 

Q.4000E-02 a { 2} 


0 . 5000E-02 a ( 3) 

-.1000E+01 end of input 
0 1 = input, 0 = stop 

P P(lst col.): to print 

0 l:to resume, 0: stop 


OUTPUT FILE: Exam3.out 

ELASTIC-PLASTIC ANALYSIS FOR CRITICAL CRACK/LOAD FOR EC02 


DATE: 17-SEP-03 TIME: 16:04:02 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
SI units [mm, MPa, MPa sqrt(mm)] 


Input Filename = Exam3.inp 
Output Filename = Exam3.out 

Thickness, t * 0.0500 

Width, W = 0.1500 

X Offset, XD = 0.0000 

[Note: Solution accurate if 2c/W < or = 0.5] 

Material Yield Stress = 400.00 

Material Ultimate Stress = 600 . 00 


Data for the 

Nonlinear Material Behavior 

SigO 

= 0.4000E+03 

E 

= 0.2100E+06 

nu 

= 0 ,3 Q0QE+Q0 

alpha 

= 0 . 1500E+01 

n 

= 0 . 1500E+02 


Data for the Elastic Plastic Failure Analysis 
* DUCTILE ANALYSIS* is performed 
Ultimate Tensile Stress (Su) = 0.6000E+03 
Jmat = 0 . 5000E-01 

Kmat(c) = 0 . 1074E+03 , Kmat(a) = 0.1074E+03 
Search for ‘CRITICAL LOAD* 

Constant aspect ratio = 0.2500E+00 

Fracture resistance curve (power law) : 

Jr = (0 . 5000E+01) *x A (0.5000E+00) 
da (max) = 0.4000E-02 

Model Code= EC02 
PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

0.00 0 . 4000E+03 

1.00 0 . 4000E+03 


TABLE OUTPUT (Fixed a/c=0 , 250E+00) : 


a P_init p_inst da (tear) dc(tear) P_init/PO 

P_inst/PO 

0 . 300E-02 0 . 316E+01 0.336E+01 0.168E-02 0.237E-03 0.107E+01 
0 . 114E+01 

0 . 400E-02 0.302E+01 0.321E+01 0.150E-02 0.167E-03 0.104E+01 
0 . 110E+01 

0 . 500E-02 0 . 288E+01 0.306E+01 0.135E-02 0.128E-03 0.100E+01 
0 . 106E+01 


EXAMPLE 4 


INPUT FILE: Exam4.inp 

Exam4 . out Output file name*12 

1 1=US units; 2=SI units 
cc Crack Model Type 

1 Crack Model Number 

1 l=tension, 2=bending 

0 . 2000E+01 Thickness 

0 . 2000E+01 Width 

0.3000E+05 Elastic Young" s modulus 

0.300 Poisson"s ratio 

1.000 Alpha 

0.1000E+03 SigmaO 

10.000 n 

0.1000E+03 material yield stress 
0.3000E+03 material ultimate stress 
2 1: p (max) , 2*.p+s(max) 

0 . 8000E+02 loading stress 

0 . 0000E+00 Non-Dimensional position 
0.5000E+02 Stress value 
0.1000E+01 Non-Dimensional position 
0.5000E+02 Stress value 
- . 1000E+01 Non-Dimensional position 
2 1: p (min) , 2:p+s(min) 

0.0000E+00 loading stress 

0.0000E+00 Non-Dimensional position 
Q.5000E+02 Stress value 
0.1000E+01 Non-Dimensional position 
0.5000E+02 Stress value 
- . 1000E+01 Non-Dimensional position 
2 no. of cycles 
2 Is p (max) , 2:p+s(max) 

0.9000E+02 loading stress 

0.0000E+00 Non-Dimensional position 
0.5000E+02 Stress value 
0.1000E+01 Non-Dimensional position 
0.5000E+02 Stress value 
- . 1000E+01 Non-Dimensional position 
2 1: p (min) , 2:p+s(min) 

0 . 00 00E+00 loading stress 


O.OOOOE+OO Non-Dimensional position 
0.5000E+02 Stress value 
0.1000E+01 Non-Dimensional position 
0.5000E+02 Stress value 
-.1000E+01 Non-Dimensional position 
1 no. of cycles 
-1 terminate input 
1 Block Case ID. 


2 no. of times 
2 Block Case ID. 

1 no. of times 


-1 Block Case ID. 


0 . 5000E+00 
0 . 1000E-08 
0 . 4000E+01 
0 . 1000E+01 
0 . 1000E+01 
0 . 1000E+01 
0 . 5000E-01 
0 . 1000E+00 
100 
10 


Jmat 

C in Paris Law 
m in Paris Law 
baseline U0 
alp_bury 
alp_surf 
initial a 
initial c 

max. no. of schedule 
print interval 


P P ( 1st col . ) : to print 
0 l:to resume, 0: stop 


OUTPUT FILE: Exam4.out 

ELASTIC -PLASTIC FATIGUE LIFE CALCULATION FOR CC01 


DATE: 18-SEP-03 TIME: 13:29:02 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
U.S . customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = exam4.inp 
Output Filename = Exam4.out 

Plate Thickness, t = 2.0000 

Plate Width, W = 2.0000 

Material Yield Stress = 100.00 

Material Ultimate Stress = 300.00 

Data for the Nonlinear Material Behavior: 


SigO 

5= 

0 . 1000E+03 

E 

5= 

0.3000E+05 

nu 

as 

0 . 3 OOOE+OO 

alpha 

as 

0 . 1000E+01 

n 

as 

0 . 1000E+02 


Model Code= CC01 under uniform tension 


BLOCK CASE DEFINITION: 

Blk Cse. Maximum Load Values Minimum Load Values, 

1 Primary Load Primary Load 

i Si i Si 

0 0 . 8000E+02 0 0.1000E-09 

1 0 . OOOOE+OO 1 0.0000E+00 

2 0. OOOOE+OO 2 0. 0000E+00 

3 0. OOOOE+OO 3 0. 0000E+00 

Secondary Load Secondary Load 

Norm, x S(Norm. x) Norm, x S(Norm. x) 

0. OOOOE+OO 0 . 5000E+02 0. OOOOE+OO 0.5000E+02 

0 . 1000E+01 0 . 5000E+02 0.1000E+01 0.5000E+02 

2 Primary Load Primary Load 

i Si i si 

0 0 . 9000E+02 0 0 . 1000E-09 

1 0. OOOOE+OO 1 0. 0000E+00 

2 0. OOOOE+OO 2 0. OOOOE+OO 

3 0. OOOOE+OO 3 0. 0000E+00 

Secondary Load Secondary Load 

Norm, x S(Norm. x) Norm, x S(Norm. x) 

0. 0000E+00 0 . 5000E+02 0. 0000E+00 0.5000E+02 

0 . 1000E+01 0 . 5000E+02 0.1000E+01 0.5000E+02 


Cycle 

2 


DEFINITION OF LOAD SPECTRUM (or SCHEDULE) : 

Blk. Case ID No. of Times Applied 

1 2 

2 1 


FATIGUE DATA (da/dN=C*dK~m) : 

J (mat) = 0 . 5000E+00 

C = 0 . 1000E-08 | 

m = 0 . 4000E+01 I 

U0 = 0 . 1000E+01 

alp (c)= 0 . 1000E+01 

alp (a)= 0 . 1000E+01 

a(init)= 0.5000E-01 

c (init) = 0 . 1000E+00 

Max. No. of Schedules = 100 

Pint interval = 10 schedules 

RESUTLS OF FATIGUE ANALYSIS: 

SCHEDULE a c D_Jeff(a) D_Jeff(c) Jmax(a) 

Jmax(c) U(a) U(c) 

10 0 . 6651E-01 0 . 1037E+00 0.2474E-01 0.1291E-01 0.1423E+00 
0 . 1074E+00 0.83 0.72 

20 0 . 8660E-01 0.1113E+00 0.2710E-01 0.1830E-01 0.1644E+00 
0 . 1392E+00 0.81 0.74 

30 0 . 1115E+00 0 . 1259E+00 0.3054E-01 0.2497E-01 0.1940E+00 
0 . 1785E+00 0.79 0.76 

40 0 . 1454E+00 0.1524E+00 0.3671E-01 0.3385E-01 0.2410E+00 
0 . 2332E+00 0.78 0.77 

50 0 . 1999E+00 0.2025E+00 0.4867E-01 0.4758E-01 0.3272E+00 
Q.3243E+QQ 0.78 0.78 


Results: two-D at c tip, J__bury (max) =0 . 5054E+00, 
J_surf (max) =0 . 5Q46E+00 

J (mat) =0 . 5000E+00 

with a=0 . 2978E+00, c=0.2985E+00 after 
total no, of cycles= 295 


58-th schedule 


EXAMPLE S 


INPUT FILE: Exam5.inp 

Exam5.out Output file name*12 

1 1=US units; 2=SI units 
1 Proof Test Procedure 

1 Proof Test - Safe Pre-Proof Life Analysis 

sc Crack Model Type 

2 Crack Model Number 

0.1000E+01 Thickness 

0 . 5000E+01 Width 

0.3000E+05 Elastic Young’s modulus 

0.300 Poisson’s ratio 

1.000 Alpha 

0.1000E+03 SigmaO 

10.000 n 

0 . 1000E+03 material yield stress 
0.3000E+03 material ultimate stress 

1 Is const, asp. 2: const, length 
0.6000E+00 constant aspect ratio 

2 1: brittle, 2: ductile 

0.2500E+00 matl toughness 

1 1: quad. 2: power 

0.2450E+00 djO — quadratic 

0.3000E+02 djl — quadratic 

- . 5000E+02 dj2 — quadratic 

0.3000E+00 da (max) — quadratic 

1 1= interactively input, 2= create a table 

2 1= pri., 2=pri. & sec. 

0 . OOOOOOOOOOOOOOOE+OOO Nondim position 

110.000000000000 Stress value 

1.00000000000000 Nondim position 

120.000000000000 Stress value 

-1.00000000000000 Nondim position 

O.OOOOE+OO Non-Dimensional position 
0.1000E+03 Stress value 
0 . 1000E+01 Non-Dimensional position 
0.5000E+02 Stress value 
-.1000E+01 Non-Dimensional position 

1.000 Load Factor # 1 

0 1 = input, 0 = stop 

P P(lst col.): to print 

0 l:to resume, 0: stop 


OUTPUT FILE: Exam5.out 

**************************************** 

* SAFE LIFE PRE-PROOF TEST ANALYSIS * 
**************************************** 


CRITICAL FLAW SIZE 


Elastic-Plastic fracture mechanics will be used to determine the 
maximum tolerable crack size that could just survive service 
conditions 


ELASTIC -PLASTIC PROOF LOAD ANALYSIS FOR SC02 


DATE: 17-SEP-03 TIME: 16:08:25 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = exam5 . inp 

Output Filename = Exam5.out 

Plate Thickness, t = 1.0000 

■ Width, W = 5.0000 

Material Yield Stress = 100.00 

Material Ultimate Stress = 300.00 

Data for the Nonlinear Material Behavior: 

SigO = 0 . 1000E+03 

E = 0 . 3000E+05 

nu = 0.3000E+00 

alpha = 0 . 1000E+01 

n = 0 . 1000E+02 

Data for the Elastic Plastic Failure Analysis 
♦DUCTILE ANALYSIS* is performed 
Ultimate Tensile Stress (Su) = 0.3000E+03 
Jmat = 0.2500E+00 

Kmat(c) = 0 . 8660E+02 , Kmat(a) = 0.9078E+02 

Search for ‘CRITICAL CRACK LENGTH* 

Constant aspect ratio = 0.6000E+00 

Fracture resistance curve (quadratic form): 

Jr = (0.2450E+00)+(0 . 3000E+02) *x+ (- . 5000E+02) *x A 2 

da(max) = 0.3000E+00 

Model Code= SC02 

PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

0.00 0 . 1100E+03 

1.00 0 . 1200E+03 

SECONDARY LOAD DISTRIBUTION: 

Norm, x Stress 

0.00 0 . 1000E+03 

1.00 0 . 5000E+02 


TABLE OUTPUT (Fixed a/c=0 . 600E+00 ) : 



* IMPORTANT NOTE * 


The Critical Flaw Size for Ductile Materials is the calculated 
critical 

crack size and equals the instability crack size minus ductile 
tearing. A 

crack with an initial size greater than the critical size will tear 
to 

instability under the given service load. 


ERROR[JSC02] : exceeds plastic collapse load! 
a=0 . 182E+00, c=0 . 304E+00 
p=0.198E+03, pmax=0 . 197E+03 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0.251E+00, c=0.418E+00 
p=0.195E+03, pmax=0 . 193E+03 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0 . 199E+00, c=0.332E+00 
p=0.198E+03, pmax=0 . 196E+03 


ERROR [JSC02] : exceeds plastic collapse load! I 
a=0 . 212E+00 , c=0 . 354E+00 ! 
p=0.197E+03, pmax=0 . 195E+03 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0 . 222E+00 , c=0.370E+00 
p=0.197E+03, pmax=0.195E+03 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0.215E+00, c=0.358E+00 
p=0.197E+03, pmax=0 . 195E+03 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0 . 216E+00 , c=0.361E+00 
p=0.197E+03, pmax=0 . 195E+03 

Pri. Load a_init a_crit a_inst da (tear) P/P0_init P/P0_crit 

P/P0_inst 

c_init c_crit c_inst dc(tear) 

0 . 100E+01 0 . 319E-01 0.216E+00 0.359E+00 0.142E+00 0.113E+01 0.116E+01 
0.120E+01 




0 . 531E-01 0 . 361E+00 0.489E+00 0.129E+00 


INPUT FILE: Exam6.inp 

Exam6.out Output file name*12 

2 1=US units; 2=SI units 
2 Proof Test Procedure 

2 Proof Test - Proof Test Analysis 

sc Crack Model Type 

4 Crack Model Number 
0 . 840000E-01 T 

0 . 118300E+01 Outer Diameter 

e 

0.2100E+06 Elastic Young's modulus 

0.300 Poisson's ratio 
1.677 Alpha 
0.5560E+03 SigmaO 

17.700 n 

0.5560E+03 material yield stress 
0.9000E+03 material ultimate stress 
0.6100E+00 constant aspect ratio 

2 Is brittle, 2: ductile 

0.4330E+00 matl toughness 

1 1: quad. 2: power 

0.3610E+00 djO — quadratic 

0.5033E+03 djl — quadratic 

-.2325E+04 dj2 — quadratic 

0.5000E-02 da (max) — quadratic 

1 1= interactively input, 2= create a table 

1 1= pri., 2=pri. & sec. 

y Interal Pressure (Y/N) 

80.450 Internal Pressure 

1.000 Load Factor # 1 

0 1 = input, 0 = stop 

P P(lst col.): to print 

0 lsto resume, 0: stop 


EXAMPLE 6 


OUTPUT FILE: Exam6.out 

**************************************** 

* PROOF TEST ANALYSIS * 

**************************************** 

FLAW SCREENING ANALYSIS 


Elastic-Plastic fracture mechanics is used to determine the 
maximum flaw 

size that could just survive the specified proof load 



ELASTIC -PLASTIC PROOF LOAD ANALYSIS FOR SC 04 


DATE: 17-SEP-03 TIME: 16:10:41 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
, SI units [mm, MPa, MPa sqrt (mm) ] 


Input Filename = exam6.inp 
Output Filename = Exam6.out 

Cylinder Thickness, t = 0.0840 

Outer Diameter, D = 1.1830 

Crack Type = EXTERNAL 

Material Yield Stress = 556.00 

Material Ultimate Stress = 900.00 

Data for the Nonlinear Material Behavior: 

SigO = 0.5560E+03 

E = 0 . 2100E+06 

nu = 0.3000E+00 

alpha = 0 . 1677E+01 

n = 0 . 1770E+02 

Data for the Elastic Plastic Failure Analysis 
‘DUCTILE ANALYSIS* is performed 
Ultimate Tensile Stress (Su) = 0.9000E+03 
Jmat = 0 . 4330E+00 I 

Kmat(c) = 0 . 3015E+03 , Kmat(a) = 0.3161E+03I 

Search for ‘CRITICAL CRACK LENGTH* 

Constant aspect ratio = 0.6100E+00 

Fracture resistance curve (quadratic form) : 

Jr = (0.3610E+00)+(0.5033E+03) *x+ (- . 2325E+04) *x A 2 
da(max) = 0.5000E-02 


Model Code= SC04 
Internal Pressure 80.450 
PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

Stresses at 10 points 


Distance x/t 
(from inner wall) 

0.0000 

0.1111 

0.2222 

0.3333 

0.4444 


Stress due to 
Int . Pressure 

529.3516 

518.4388 

508.1016 

498.3002 

488.9981 


Total Stress 


529.3516 
518.4388 
508.1016 
498.3002 
488 . 9981 



0.5556 

480.1621 

480.1621 

0.6667 

471.7615 

471.7615 

0.7778 

463.7681 

463.7681 

0.8889 

456.1561 

456.1561 

1.0000 

448.9016 

448.9016 


NO SECONDARY LOAD SPECIFIED! 


TABLE OUTPUT (Fixed a/c=0 . 610E+00) : 


* IMPORTANT NOTE * 


The flaw size for Ductile Materials that will just survive the Proof 
Load is the calculated instability crack size and equals the crtical 
crack size plus ductile tearing. A crack with an initial size 
greater 

than the size to initial ductile tearing but less than the critical 
size will tear under the Proof Load but not fail 


Pri . Load a_init a_crit 

P/P0_inst 

c_init c_crit 

0 . 100E+01 0 . 562E-01 0.610E-01 
0 . 781E+01 


a_inst da (tear) 

c_inst dc(tear) 

0 . 634E-01 0 . 241E-02 


P/P0_init P/P0_crit 
0 . 723E+01 0 . 757E+01 


0 . 921E-01 0 . 100E+00 


0 . 104E+00 0 . 352E-02 


EXAMPLE 7 


INPUT FILE: Exam7.inp 

ex7.out Output file name*12 

1 1=US units; 2=SI units 

2 Proof Test Procedure 

1 Proof Test - Proof Test Analysis 

ec Crack Model Type 

2 Crack Model Number 

0 . 500E+00 Thickness 
0 . 600E+01 Width 

0.3000E+05 Elastic Young's modulus 

0.300 Poisson's ratio 

1.000 Alpha 

0.8000E+02 SigmaO 

25.000 n 

0 . 8000E+02 material yield stress 
0 . 1200E+03 material ultimate stress 


0.2500E+00 constant aspect ratio 
2 Is brittle, 2: ductile 

0.2000E+00 matl toughness 

2 1: quad. 2: power 

0.5000E+01 djl -- power law 

0 . 5000E+00 dj2 — power law 

0.10000E+00 da (max) — power law 
1 1: with 2nd load, 0: w/o 2nd load 


O.OOOOOOOOOOQOOOOE+OOO 
60.0000000000000 
0.200000000000000 

55.0000000000000 
0.400000000000000 

52.0000000000000 
0.600000000000000 

50.0000000000000 
0.800000000000000 

49.0000000000000 

1.00000000000000 
49.0000000000000 

- 1.00000000000000 
0.0000E+00 Non-Dimensional position 
0.9000E+02 Stress value 
0.1000E+01 Non-Dimensional position 
O.OOOOE+OO Stress value 
-.1000E+01 Non-Dimensional position 
1 1= interactively input, 2=create a table 

0 . 2500E-01 a( 1) 

0 . 3000E-01 a( 2) 


Nondim position 
Stress value 
Nondim position 
Stress value 
Nondim position 
Stress value 
Nondim position 
Stress value 
Nondim position 
Stress value 
Nondim position 
Stress value 
Nondim position 


0 . 3500E-01 
0 . 4000E-01 
0 . 4500E-01 
0 . 5000E-01 
- . 1000E+01 
P P(lst 
0 1 : to 


a< 3) 
a ( 4) 
a ( 5) 
a( 6) 

end of input 
col . ) s to print 
resume, 0: stop 


l 


OUTPUT FILE: Exam7.out 

**************************************** 

* PROOF TEST ANALYSIS * 

**************************************** 

PROOF LOAD ANALYSIS 


Elastic-Plastic fracture mechanics will be used to determine the 
proof load necessary to screen against the presence of specified 
initial crack sizes 


ELASTIC - PLASTIC PROOF LOAD ANALYSIS FOR EC02 




DATE: 18-SEP-03 TIME: 22:58:27 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = ex7.inp 
Output Filename = ex7.out 

Thickness, t = 0.5000 

Width, W = 6.0000 

X Offset, XD = 0.0000 

[Note: Solution accurate if 2c/W < or = 0.5] 

Material Yield Stress = 80.00 

Material Ultimate Stress = 120.00 


Data for the 

Nonlinear Material Behavior: 

SigO 

= 0 . 8000E+02 

E 

= 0 . 3000E+05 

nu 

= 0 . 3000E+00 

alpha 

= 0 . 1000E+01 

n 

= 0 . 2500E+02 


Data for the Elastic Plastic Failure Analysis 
♦DUCTILE ANALYSIS* is performed 
Ultimate Tensile Stress (Su) = 0.1200E+03 
Jmat = 0.2000E+00 

Kmat(c) = 0 . 8120E+02 , Kmat(a) = 0.8120E+02 
Search for ‘CRITICAL LOAD* 

Constant aspect ratio = 0.2500E+00 

Fracture resistance curve (power law) : 

Jr = (0 . 5000E+01) *x A (0.5000E+00) 
da (max) = 0.1000E+00 


Model Code= EC02 
PRIMARY LOAD DISTRIBUTION 1: 

Norm, x 
0.00 
0.20 
0.40 
0.60 
0.80 
1.00 

SECONDARY LOAD DISTRIBUTION: 

Norm, x Stress 

0.00 0 . 9000E+02 

1.00 0 . OOOOE+OO 


Stress 
0 . 6000E+02 
0.5500E+02 
0.5200E+02 
0 . 5000E+02 
0.4900E+02 
0 . 4900E+02 


* IMPORTANT NOTE * 






The Proof Load Factor for a Ductile Material containing a given 
initial 

crack size is equal to the instability load factor. A proof load 
that 

exceeds the load to initiate tearing but is less than the 
instability 

load will cause a crack to tear but not fail. 


TABLE OUTPUT (Fixed a/c=0 . 250E+00) : 

a P_init P_inst da (tear) dc(tear) P_init/P0 

P_inst/P0 

0 . 250E-01 0 . 260E+03 0.270E+03 0.353E-01 0.943E-02 0.109E+01 
0 . 113E+01 

0 . 300E-01 0 . 257E+03 0.267E+03 0.349E-01 0.791E-02 0.107E+01 
0 . 112E+01 

0 . 350E-01 0 . 254E+03 0.265E+03 0.340E-01 0.659E-02 0.106E+01 
0 . 111E+01 

0 . 400E-01 0 . 251E+03 0.262E+03 0.328E-01 0.555E-02 0.105E+01 
0.110E+01 

0 . 450E-01 0 . 249E+03 0.260E+03 0.317E-01 0.473E-02 0.104E+01 
0 . 109E+01 

0 . 500E-01 0 . 246E+03 0.258E+03 0.306E-01 0.410E-02 0.104E+01 
0 . 108E+01 

EXAMPLE81 


INPUT FILE: Exam8.inp 

Exam8.out Output file name*12 

1 1=US units; 2=SI units 

2 Proof Test Procedure 

3 Proof Test - Proof Test Analysis 

cc Crack Model Type 

1 Crack Model Number 

2 l=tension, 2=bending 

0 . 2000E+01 Thickness 

0 . 2000E+01 Width 

0.3000E+05 Elastic Young’s modulus 

0.300 Poisson’s ratio 
1.200 Alpha 
0.1500E+03 SigmaO 

20.000 n 

0 . 1500E+03 material yield stress 
0.2000E+03 material ultimate stress 
0 . 5000E+00 Initial Proof Crack Length 

0.1000E+01 Init Proof a/c 

0 proof secondary (0-no 1-yes) 

0 . 2200E+03 loading stress 



1 1: quad. 

0 . 1450E+00 
0 . 3000E+02 
- . 5000E+02 
0 . 3000E+00 
0 . 1500E+00 


2 : power 

djO — quadratic 
djl — quadratic 
dj2 — quadratic 
da (max) — quadratic 
Jmat 


P P(lst col.): to print 


0 1 : to resume, 0: stop 


OUTPUT FILE: Exam8.out 

*******★*******★*★******★★***★*★★***★★** 

* PROOF TEST ANALYSIS * 

**************************************** 

FINAL FLAW SIZE ANALYSIS 


Elastic-Plastic fracture mechanics is used to determine the final 

flaw 

size after the application of the specified proof load 


ELASTIC -PLASTIC PROOF LOAD ANALYSIS FOR CC01 


DATE: 17-SEP-03 TIME: 16:14:11 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Aug. 2002 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = exam8 . inp 
Output Filename = Exam8.out 

Plate Thickness, t = 2.0000 

Plate Width, W = 2.0000 

Material Yield Stress = 150.00 

Material Ultimate Stress = 200.00 


Data for the 

Nonlinear Material tsehavior: 

SigO 

= 0 . 1500E+03 

E 

= 0 . 3000E+05 

nu 

= 0 . 3000E+00 

alpha 

= 0 . 1200E+01 

n 

= 0 . 2000E+02 


PRIMARY LOAD DISTRIBUTION 

0.0000 SO: Tensile Stress 


220.0000 


SI: Bending Stress 



0.0000 


S2 : Bending Stress 


Fracture resistance curve (quadratic form): 

Jr = (0 .1450E+00)+(0.3000E+02) *x+ (- . 5000E+02) *x A 2 
da (max) = 0.3000E+00 

, Toughness =0.1500E+00 

END OF PROOF TEST - FINAL FLAW SIZE ANALYSIS: 


a 

c 

da 

dc 

0 . 5694E+00 

0 . 6555E+00 

0.6947E-01 

0 . 1557E+00 



EXAMPLE 9 


INPUT FILE: Exam9.inp 

Exam9.out Output file name*12 

1 1=US units; 2=SI units 
sc Crack Model Type 

4 Crack Model Number 

0 . 500000E+00 T 

0 . 400000E+02 Outer Diameter 

e 

0 . 3000E+05 Elastic Young’s modulus 

0.300 Poisson’s ratio 

1.000 Alpha 

0 . 1000E+03 SigmaO 

10.000 n 

0.1000E+03 material yield stress 
0 . 1500E+03 material ultimate stress 
0 . 2500E+00 Initial Proof Crack Length 

0 . 5000E+00 Init Proof a/c 

100 Number of Service Cycles 

1 secondary (0-no 1-yes) 
n 

0 . OOOOOOOOOOOOOOOE+OOO Nondim position 

90.0000000000000 Stress value 

1.00000000000000 Nondim position 

80.0000000000000 Stress value 

-1.00000000000000 Nondim position 

0.0000E+00 Non-Dimensional position 
-0.100E+03 Stress value 
0.1000E+01 Non-Dimensional position 
0.1000E+03 Stress value 
-.1000E+01 Non-Dimensional position 
n 

0. OOOOOOOOOOOOOOOE+OOO Nondim position 
0 . OOOOOOOOOOOOOOOE+OOO Stress value 

1.00000000000000 Nondim position 

0 . OOOOOOOOOOOOOOOE+OOO Stress value 

-1.00000000000000 Nondim position 

0.0000E+00 Non-Dimensional position 
-0.100E+03 Stress value 
0.1000E+01 Non-Dimensional position 
0.1000E+03 Stress value 
-.1000E+01 Non-Dimensional position 
1 1 : quad . 2 : power 

0.2000E+00 djO — quadratic 

0 . 3000E+02 djl — quadratic 

0 . 0000E+00 dj2 — quadratic 

0.2000E+00 da (max) — quadratic 

0 . 2100E+00 Jmat 

0 . 1000E+01 baseline U0 

0 . 1000E+01 alp_bury 

0.1000E+01 aip_surf 



0.1000E-09 C in Paris Law (Serv) 

0.4000E+01 m in Paris Law (Serv) 

P P(lst col.): to print 

0 l:to resume, 0: stop 


OUTPUT FILE: Exam9.out 


TEAR-FATIGUE ANALYSIS FOR SC04 


DATE: 19-SEP-03 TIME: 13:53:44 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Jan. 2003 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = exam9.inp 
Output Filename = Exam9.out 

Cylinder Thickness, t = 0.5000 

Outer Diameter, D = 40.0000 

Crack Type = EXTERNAL 

Material Yield Stress = 100.00 

Material Ultimate Stress = 150.00 

Data for the Nonlinear Material Behavior: 
SigO = 0 . 1000E+03 
E = 0 . 3000E+05 

nu = 0.3000E+00 

alpha = 0.1000E+01 
n = 0 . 1000E+02 


SERVICE SPECTRUM 


LOAD STEP 1 


CYCLIC MAXIMUM DISTRIBUTION 


PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

Stresses at 2 points 


Distance x/t Stress due to 

(from inner wall) Int. Pressure 


0.00 0 . 9000E+02 

1.00 0 . 8000E+02 

SECONDARY LOAD DISTRIBUTION: 

Norm, x Stress 


Total Stress 



0.00 

1.00 


0 . 1000E+03 
- . 1000E+03 


CYCLIC MINIMUM DISTRIBUTION 


PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

Stresses at 2 points 


Distance x/t Stress due to 

(from inner wall) Int. Pressure 


0.00 0 . OOOOE+OO 

1.00 0. 0000E+00 

SECONDARY LOAD DISTRIBUTION: 

Norm, x Stress 

0.00 0 . 1000E+03 

1.00 - . 1000E+03 


Total Stress 


Fracture resistance curve (quadratic form): 

Jr = ( 0 . 2000E+00) + (0 . 3000E+02) *x+ (0 . 0000E+00) *x A 2 
da (max) = 0.2000E+00 


FATIGUE DATA (da/dN=C*dK~m) : 

J (mat) = 0 . 2100E+00 
C = 0 . 1000E-09 

m = 0 . 4000E+01 

U0 = 0 . 1000E+01 

alp (c)= 0 . 1000E+01 
alp (a) = 0 . 1000E+01 

a(init)= 0.2500E+00 
c(init)= 0.5000E+00 

RESULTS OF TEAR-FATIGUE ANALYSIS: 


Cycle 

a 

c 

da 

dc 

0 

0 . 2500E+00 

0 . 5000E+00 



0 

0 . 2578E+00 

0 . 5133E+00 

0.7968E-02 

0 . 1365E-01 

1 

0 . 2729E+00 

0 . 5310E+00 

0 . 1509E-01 

0 . 1767E-01 

2 

0 . 2879E+00 

0 . 5514E+00 

0 . 1502E-01 

0.2045E-01 

3 

0.3029E+00 

0 . 5750E+00 

0 . 1499E-01 

0.2356E-01 

4 

0.3179E+00 

0 . 6021E+00 

0 . 1502E-01 

0.2706E-01 

5 

0 . 3330E+00 

0 . 6331E+00 

0 . 1513E-01 

0 . 3102E-01 

6 

0 . 3484E+00 

0 . 6686E+00 

0 . 1533E-01 

0.3557E-01 

7 

0 . 3640E+00 

0.7095E+00 

0 . 1566E-01 

0 . 4086E-01 

8 

0 . 3802E+00 

0 . 7566E+00 

0.1615E-01 

0.4713E-01 

9 

0 . 3971E+00 

0 . 8114E+00 

0 . 1693E-01 

0 . 5481E-01 

10 

0.4156E+00 

0 . 8764E+00 

0 . 1845E-01 

0 . 6493E-01 

11 

0 . 4390E+00 

0 . 9585E+00 

0 . 2344E-01 

0 . 8210E-01 

Failed 

on cycle 

12 




END OF TEAR- FATIGUE ANALYSIS: 




Cycle 


a 


c 


11 0 . 43 90E+00 0 . 9585E+00 


EXAMPLE 10 


MULT I -CYCLE PROOF TEST ANALYSIS FOR SCO 2 


DATE: 23 -SEP-03 TIME: 15:51:23 
(computed: NASA/FLAGRO Version 3.00, October 1995.) 
Elastic-Plastic Fracture Module (EPFM) V.x.xx, Jan. 2003 
U.S. customary units [inches, ksi, ksi sqrt(in)] 


Input Filename = examlO.inp 
Output Filename - ExamlO.out 

Plate Thickness, t = 1.0000 

■ Width, W « 10.0000 

Material Yield Stress = 100.00 

Material Ultimate Stress « 200.00 

Data for the Nonlinear Material Behavior: 


SigO 

K 

0.1000E+03 

E 

= 

0 . 3 000E+05 

nu 

= 

0 . 3 OOOE+OO 

alpha 

= 

0 . 1000E+Q1 

n 


0 . 5000E+01 

PROOF LOAD 




Number of Proof Cycles 


4 


CYCLIC MAXIMUM DISTRIBUTION 

PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

0.00 0 . 1200E+03 

1.00 0 . 1200E+03 

CYCLIC MINIMUM DISTRIBUTION 


PRIMARY LOAD DISTRIBUTION 1: 


Norm, x 
0.00 
1.00 


Stress 
0 . 0000E+00 
0 . 0000E+00 


SERVICE SPECTRUM 


Number of Service Cycles 


50 


LOAD STEP 1 

Load Step 1 applied 2 cycles per service spectrum 


CYCLIC MAXIMUM DISTRIBUTION 

PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

0.00 0 . 8000E+02 

1.00 0.8000E+02 

CYCLIC MINIMUM DISTRIBUTION 

PRIMARY LOAD DISTRIBUTION 1: 

Norm, x Stress 

0.00 0 . 0000E+00 

1.00 0.0000E+00 

Fracture resistance curve (quadratic form) : 

Jr = (0 . 2000E+00) + (0 . 2000E+02) *x+ (0 . 0000E+00) *x A 2 
da (max) = 0.2000E+00 


PROOF CYCLE FATIGUE DATA 


FATIGUE DATA (da/dN=C*dK A m) : 
J(mat) = 0.2100E+00 
C = 0 . 1000E-08 

m = 0 . 4000E+01 

U0 = 0 . 1000E+01 

alp (c) = 0 . 1000E+01 
alp (a) = 0 . 1000E+01 

a(init)= 0.0000E+00 
c(init)= 0 . OOOOE+OO 

SERVICE CYCLE FATIGUE DATA 


FATIGUE DATA (da/dN=C*dK A m) : 
J (mat) = 0 . 2100E+00 
C = 0.2000E-08 

m = 0 . 4000E+01 



UO = 0 . 1000E+01 

alp (c)= 0 . 1000E+01 
alp (a)« 0 . 1000E+01 

a(init)= O.OOOOE+OO 
c(init)= O.OOOOE+OO 

The MCPT module employs a root -finding algorithm to determine the 
initial 

crack lengths: 

Proof Load Only : maximum initial crack length which survives 
prescribed 

number of Proof Load applications. 

Service Load : minimum initial crack length which survives 
prescribed 

number of Proof Load and Service Load applications. 


During the initial crack length searches various error and warning 
messages 

may be displayed. These messages do not affect the probability 
calculation. 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0 . 103E+01 , c=0 . 133E+01 
p=0 . 120E+03 , pmax=0 . 118E+03 


ERROR [JSC02] : exceeds plastic collapse load! 
a=0 . 102E+01, c=0 . 131E+01 
p=0 . 120E+03 , pmax=0 . 119E+03 


MCPT ANALYSIS COMPLETE 


Mean of Exponential Dist 0.150000000000000 

Conditional Probability of Failure 0.254858919567750 

Probability of Proof Failure 0.806209301940743 

Probability of Proof and Service Failure... 0.600739670302655 

Initial crack size to just survive proof test 

0.246146486746147 

Initial crack size to just survive proof test and service . . . . 
0.137721242920040 



APPENDIX 4: VALIDATION AND CONSISTENCY CHECKS 


This appendix lists in tabular form the results of part of the exercise performed to validate 
the NASGRO EPFM and Proof Test Analysis Modules (hereafter referred to as the Modules). 
Except for the Option 5 (J estimation), all the analyses used in the validation involved 2-DOF, 
and except for Option 5 and Option 7 (fatigue lifetime), all the validation analyses addressed 
ductile fracture behavior. 


Table A4.1: Validation of Example 1: J D Solutions for SC04: Options 


Ri 


2.75 


2.75 


2.75 


2.75 


2.75 


2.75 


0.25 


0.25 


0.25 


0.25 


0.25 


0.25 


10 


10 


10 


10 


10 


10 


0.05 


0.05 


0.05 


0.1 


0.1 


0.15 


0.05 


0.1 


0.15 


0.1 


0.15 


0.15 


M 


1 .000365 


1.00146 


1.003282 


1.002917 


1.006552 


1.009813 


1 


1 


1 


1 


1 


1 


Rm 


2.875 


2.875 


2.875 


2.875 


2.875 


2.875 


yield 


100 


100 


100 


100 


100 


100 


10 


10 


10 


10 


10 


10 


mu 


0.82416 


0.82416 


0.82416 


0.82416 


0.82416 


0.82416 


Po 


9.09008 


9.087597 


9.083481 


9.073313 


9.051626 


8.960303 


Je(a) 


0.0304 


0.0563 


0.0695 


0.0621 


0.0957 


0.0954 


Je(c) 


0.0431 


0.0399 


0.0347 


0.0931 


0.0987 


0.156 


Jp(a) (manual) 


0.0591 


0.1098 


0.1361 


0.1228 


0.1934 


0.2112 


Jp(a) (NASGRO SC04) 


0.0591 


0.1100 


0.1360 


0.1230 


0.1930 


0.2110 


Jp(c) (manual) 


0.1017 0.0944 


0.0824 


0.2234 


0.2420 


0.4190 


Jp(c) (NASGRO SC04) 


0.1020 0.0943 


0.0825 


0.2230 


0.2420 


0.4180 











































































































Table A4.2: Validation: J p Solutions for SC02; Option 5 (Force, P) 


W 


alpha 


n 


mu 


Po 


Je(a, primary onl 


Je(c, primary onl 


1 

1 

320 

320 

0.2 

Km 

IHSEIB 

0.4 

1.8164 

1.8164 

1.256 

1.256 

80 

80 

1.5 

1.5 

15 

15 

0.82416 

0.82416 

314.970 

309.95 

0.0685 

0.1360 

0.1420 

0.1320 


320 


0.2 


0.8 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


299.89 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


299.89 


0.1070 0.2780 


.8 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


279.79 


0.2290 


0.3187 


320 


.8 


0.8 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


239.58 


0.2050 


0.7180 


Jp(a) (manual 


Jp(a) (NASGRO SC04 


0.192 


.192 


.478 


0.476 


1.08 

0.540 

3.37 

26.5 

1.08 

0.540 

3.37 

26.5 


Jp(c) (manual 


Jp(c) (NASGRO SC04 


0.334 


.333 


.389 

0.500 

1.30 

3.93 

77.8 

.388 

0.502 

1.30 

3.92 

77.8 






























































































Table A4.3: Validation: J« Solutions for SC02: Option 5 (Force, P, and Moment, M) 


W 


t 


Force, P, and moment, M corresponding 
to the primary linear stress 120-120(x/t) 
=60+60(1 -2(x/t 


A=M/Pt 


m 


alpha 


n 


Mu 


P=240 

M=40 


0.1667 


0.2 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


1 


314.970 


0.0868 


0.1520 


P=240 

M=40 


0.1667 


0.4 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


1 


309.95 


0.1680 


0.1410 


P=240 

M=40 


0.1667 


0.2 


0.8 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


1 


299.89 


0.2400 


0.1140 


P=240 

M=40 


0.1667 


0.4 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


1 


299.89 


0.1300 


0.3130 


P=240 

M=40 


0.1667 


0.8 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


1 


279.79 


0.3010 


0.3560 


P=240 

M=40 


0.1667 


0.8 


0.8 


1.8164 


1.256 


80 


1.5 


15 


0.82416 


1 


239.58 


0.1120 


0.8210 


Jp(a) (manual 


Jp(a) (NASGRO SC04 


0.430 

1.11 

2.89 

1.488 

11.7 

41.0 

0.430 

1.11 

2.89 

1.470 

11.7 

40.8 


Jp(c) (manual 


Jp(c) (NASGRO SC04 


0.646 

0.799 

1.18 

3.06 

11.9 

258 

0.645 

0.799 

1.17 

3.05 

11.9 

258 


















































































































Table A4.4: Validation: J D 

Solutions for CC01 Option 5 (Moment, M) 


W 

2 

2 

2 

2 

2 

2 

t 

2 

2 

2 

2 

2 

2 

Moment, M, corresponding to primary bending 
stress 1 60 

320 

320 

320 

320 

320 

320 


alpha 


n 


mu 


Mo 


0.2 


0.8089 


0.8089 


80 


1.5 


15 


1 


0.8089 


0.8089 


80 


1.5 


15 


1 


314.970 309.95 


0.8 


0.8089 


0.8089 


80 


1.5 


15 


1 


299.89 


0.8089 


0.8089 


80 


1.5 


15 


1 


299.89 


0.8 


0.8089 


0.8089 


80 


1.5 


15 


1 


279.79 


0.8 


0.8 


0.8089 


0.8089 


80 


1.5 


15 


1 


239.58 


0.127 

0.208 

0.269 

0.187 

0.353 

0.223 

0.160 

0.147 

0.115 

0.310 

0.302 

0.667 


Jp(a) (manual 


Jp(a) (NASGRO SC04 


0.167 

0.340 

0.683 

0.358 

1.631 

3.07 

0.167 

0.339 

0.682 

0.357 

1.63 

3.07 


Jp(c) (manual 


Jp(c) (NASGRO SC04 


0.211 

0.240 

0.292 1 

0.593 

1.40 

9.19 

0.210 

0.240 

0.291 

0.593 

1.40 

9.19 



































































































Table A4.5: Validation: J D Solutions for EC02: OptionS 


W 


t 2 I 2 | 2 I 2 


Force, P, corresponding to primary quadratic 
stress 120-240(x/t)+240(x/tr 
Secondary stress 60-60(x/t) 320 



alpha 


n 


mu 


Po 


Jef+a, primary onl 


Je(c, primary onl 


Je(-a, primary onl 


Jp(+a) (manual 


Jp(+a) (NASGRO SC04 



2 


1.6575 


1.6575 


80 


1.5 


15 


0.82416 


625.81 


0.0307 


0.0280 


0.0307 


1.6575 


1.6575 


80 


1.5 


15 


0.82416 


600.94 


0.0507 


0.0237 


0.0507 


0.0861 

0.251 

0.0861 

0.251 


0.8 


1.6575 


1.6575 


80 


1.5 


15 


0.82416 


554.78 


0.0648 


0.0153 


0.0648 


.982 


0.981 



1.6575 


1.6575 


80 


1.5 


15 


0.82416 


565.20 


0.0719 


0.0576 


0.0719 


.8 


1.6575 


1.6575 


80 


1.5 


15 


0.82416 


482.91 


0.1200 


0.0496 


0.1200 


320 


.8 


0.8 


1 .6575 


1.6575 


80 


1.5 


15 


0.82416 


371.33 


0.2820 


0.1300 


0.2820 


0.839 

12.7 

1179 

0.839 

12.7 

1180 



Jp(c) (manual) 

0.0785 

0.117 

0.232 

0.672 

Jp(c) (NASGRO SC04) 

0.0786 

0.117 

0.232 

0.672 

Jp(-a) (manual) 
Jp(-a) (NASGRO SC04) 

0.0861 

0.0861 

0.251 

0.251 

0.982 

0.981 

0.839 

0.839 


5.24 

544 

5.24 

543 

12.7 

12.7 

1179 

1180 













































































































Table A4.6: Validation of Example 2: Critical Crack Failure Analysis 
Instability Predictions for 2-DOF Cracks for SC02: Option 6 


Resistance Curve: 0.245+30 Aa - 50 (Aa) A 2 


Instability Conditions: 

J(a) = JR(Aa) and J(c) = JR(Ac) 
X-A)(Y-D) - BC = 0 


Derivatives 


0.3342 

0.641 

3.56 

2.56 

0.341 

0.641 

3.62 

2.66 

0.3478 

0.641 

3.69 

2.76 

0.3342 

0.6282 

3.53 

2.63 

0.3342 

0.641 

3.62 

2.66 

0.3342 

0.6538 

3.71 

2.68 


7.020 


14.663 


1.950 


0.151 0.096 


The accuracy of the derivatives 


is limited by the fact that NASGRO 


outputs are only printed to three 


significant figures. 




99.06 


BC 


102.94 




Table A4.7: Validation of Example 2: Consistency Check: SC02 


Critical Crack Failure Analysis: Option 6 


a- 


0.191 


0.545 


Initial crack 
size, a-tip 

Initial crack 
size, c-tip 

0.191 

0.545 


Applied primary 
load 

Applied secondary 
load 

1 

iio 

1 100-50(x/t) 


Check: Critical Load Failure Analysis: Option 6 


P 


0.151 


0.0961 


Predicted 
instability 
primary load 

Applied secondary 
load 

114 

100-50(x/t) 


























































Table A4.8: Validation of Example 3: Consistency Check: EC02 


Critical Load Failure Analysis: Option 6 


Initial crack 
size, a-tip 

Initial 

aspect ratio, 
a/c 

Predicted 
tear at 
instability , 
a-tip 

Predicted 
tear at 
instability, 
c-tip 

Predicted 
instability 
primary load 

Applied secondary 
load 

0.003 

0.25 

0.00168 

0.000237 

448 

0 

0.004 

0.25 

0.0015 

0.000167 

428 

0 

0.005 

0.25 

0.00135 

0.000128 

408 

0 


Check: Critical Crack Failure Analysis: Option 6 












































Table A4.9: Validation of Example 4: 

: Consistency Check: CC01 (tension) 

Fatigue Life Analysis: Option 7 

Check: Safe Life Pre-Proof Test Fatigue 
Life Analysis: Option 8 


Schedule 


Schedule 


0 


10 


20 


30 


40 


50 


58 (Failure) 


0.05 


.06651 


.08660 


.1115 


0.1454 


0.1999 


.2978 


c- 


0.1 


0.1037 


0.1524 


0 


10 


20 


30 


40 


50 


58 (Failure) 


a- 


0.05 


0.06651 


0.08660 


0.1115 


0.1454 


0.1999 


0.2978 


0.1 


0.1037 


0.1113 


0.1259 


0.1524 


.2025 


.2985 


0.215 


Table A4.10: Validation of Example 5: Consistency Check: SC02 


Safe Life Pre-Proof Test Critical Crack Analysis: Option 8 


Predicted 
Critical 
crack size, 

c-tip a- 

0.359 0.142 I 0.129 


Check: Critical Load Failure Analysis: Option 6 


P 


Applied primary 
load 

Applied secondary 
load 

110+10(x/t) 

100-50(x/t) 

(tensile force = 

1 

575) 

i 


Initial crack 
size, a-tip 

Initial crack 
size, c-tip 

0.215 

0.359 


0.142 


0.129 


Predicted 

Applied secondary 

instability 

load 

primary load 


575 

(tensile force) 

100-50(x/t) 


























































Table A4.ll; Validation of Example 6 : Consistency Check: SC04 


Proof 

Test Analysis: Flaw Screening: Option 8 

Predicted 

Maximum 

Flaw Size 

Screened 

(crack size 

at incipient 

instability, 

includes 

tearing), 

a-tip 

Predicted 

Maximum 

Flaw Size 

Screened 

(crack size 

at incipient 

instability 

includes 

tearing), 

c-tip 

Predicted 
tear at end 
of proof 
test, 
a-tip 

Predicted 
tear at end 
of proof 
test, 
c-tip 

Applied proof 
primary load 

Secondary load 
present during 
proof test 

0.0634 

0.104 

0.00241 

0.00352 

80.45 (pressure) 

"o 

Check: Proof Test Analysis: Proof Load: Option 8 

Initial 
crack size, 
a-tip (no 
tearing 
included) 

Initial 

aspect ratio, 
a/c 

Predicted 
tear at 
instability, 
a-tip 

Predicted 
tear at 
instability, 
c-tip 

Predicted Proof 
Primary Load 
Needed to Screen 
Initial Flaw 

Applied secondary 
load 

0.061 

0.61 

0.00241 

0.00352 

80.5 (pressure) 

0 

Check: Critical Load Failure Analysis: Option 6 

Initial crack 
size, a-tip 
(no tearing 
included) 

Initial crack 
size, c-tip 
(no tearing 
included) 

Predicted 
tear at 
instability, 
a-tip 

Predicted 
tear at 
instability, 
c-tip 

Predicted 
instability 
primary load 

Applied secondary 
load 

0.061 

0.100 

0.00241 

0.00352 

80.5 (pressure) 

0 

Critical Crack Failure Analysis: Option 6 

Predicted 
Critical 
crack size, 
a-tip 

(no tearing 
included) 

Initial 

aspect ratio, 
a/c 

Predicted 
tear at 
instability, 
a-tip 

Predicted 
tear at 
instability, 
c-tip 

Applied primary 
load 

Applied secondary 
load 

0.061 

0.61 

0.00241 

0.00352 

80.45 (pressure) 

0 


This example is based on a full-scale pressure vessel test. The material properties are those 
measured on the pressure vessel steel. In the test, an axial crack was machined into the outer 
surface of the vessel with a depth of 61mm and a total surface length of 200mm. This crack 
extended by ductile tearing under increasing hydraulic pressure until the crack failed (caused a 


leak) at a pressure of 80.45 MPa. 







































Table A4.12: Validation of Example 7: Consistency Check: EC02 


Proof Test Analysis: Proof Load: Option 8 


Predicted proof 


Initial crack 
size, a-tip 

Initial 

aspect ratio, 
a/c 

0.025 

0.25 

0.035 

0.25 

0.045 

0.25 


0.0250 


0.0353 


0.0453 



primary load 


Applied secondary 
load 



0.0353 I 0.00943 


0.0340 0.00659 


0.0316 0.00473 


Check: Proof Test Analysis: Flaw Screening: Option 8 


Initial 

aspect ratio, 
a/c 


90-90(x/t) 


90-90(x/t) 




0.0353 


0.0339 


0.0316 


0.00944 


0.00652 


0.00469 


Applied primary 
load 

Applied secondary 
load 

270 

90-90(x/t) 

265 

90-90(x/t) 

260 

90-90(x/t) 


Table A4.13: Validation of Example 8: Consistency Check: CC01 (bend) 


Proof Test Analysis: 

Final Crack Size: Option 8 


Initial crack 
size, a-tip 

Initial 

aspect ratio, 
a/c 

Predicted 
crack size 
after proof 
load 
applied, 
a-tip 

Predicted 
crack size 
after proof 
load 
applied, 
c-tip 

Proof primary 
load (bend) 

Applied secondary 
load 

0.5 

0.5 

0.5694 

0.6555 

220 

0 

Check: Tear-Fatigue Analysis: Option 9 
(predicted crack extension on first load application) 

Initial crack 
size, a-tip 

Initial 

aspect ratio, 
a/c 

Predicted 
crack size 
after service 
load 
applied, 
a-tip 

Predicted 

crack size 

after 

service 

load 

applied, 

c-tip 

Service primary 
load (bend): 
maximum, 
minimum 

Applied secondary 
load 

0.5 

0.5 

0.5694 

0.6555 

220,0 

0 




























































Table A4.14: Validation of Example 9: Consistency Check: SC04 


Tear-Fatigue Analysis: Option 9 


Initial crack 
size, a-tip 

Initial 

aspect ratio, 
a/c 

Predicted 
crack size 
after service 
load 
applied, 
a-tip 

Predicted 

crack size 

after 

service 

load 

applied, 

c-tip 

Service primary 
load: maximum, 
minimum 
(stress origin 
with respect to 
inner surface) 

0.25 

0.5 

0.2578 

0.5133 

90-10(x/t), 0 


Applied secondary 
load 


100-200(x/t) 


Check: Proof Test Analysis: Final Crack Size: Option 8 


Initial crack 
size, a-tip 

Initial 

aspect ratio, 
a/c 

0.25 

0.5 


0.2578 


Predicted 
crack size 
after proof 
load 
applied, 
c-tip 

Proof primary 
load 

(stress origin 
with respect to 
inner surface) 

Applied secondary 
load 

0.5133 

90-10(x/t) 

100-200(x/t) 


Table A4.15: Validation of Example 10: SC02 


Multiple Cycle Proof Test: Option 10 


Crack size, a p , 

Crack size, a,*. 

Probability of 

Probability of 

Conditional 

to just survive 

to just survive 

failing during 

failing during the 

probability of 

the proof test 

the proof test 

the proof test: 

proof test and 

failure for 

predicted by 

and service 


service: 

specified service 

NASGRO 

predicted by 



lifetime given 

Option 10 

NASGRO 



survival of MCPT 

Analysis 

Option 10 





Analysis 

p =i-<A ois J 

p„= i-^ 0l5 > 


0.26104 

0.13763 

r 


“cond p 

r p 

Manual 

0.8245 

0.6005 

0.2717 

NASGRO SC02 

0.8245 

0.6005 

0.2717 




















Table A4.16: Validation: SC02 

Multiple Cycle Proof Test: Option 10: 
(No Proof Test) 

Check: Fatigue Life Analysis: Option 7 

Predicted Crack size, a-tip, to just survive 
service lifetime of 50 cycles 

Predicted service lifetime (cycles) for initial 
crack size of 0.21214 

0.21214 

49 


Table A4.17: Validation: SC02 

Multiple Cycle Proof Test: Option 10: 
(No Service) 

Check: Tear-Fatigue Analysis: Option 9 

Predicted Crack size, a-tip, to just survive 4 
proof cycles 

Predicted lifetime (cycles) for initial crack size 
of 0.26104 


0.26104 
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